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1.0  INTRODUCTION 


Aircraft  structural  systems  must  meet  stringent  strength  demands,  as  well  as 
requirements,  for  durability  and  integrity  during  their  operational  lifetime.  As  part  of 
the  related  research  and  development  work  for  materials,  the  Air  Force  Materials 
Laboratory  (AFML)  has  sponsored  considerable  research  concerned  with  the  application 
of  reliability  technology  to  the  problem  of  assessing  and  assuring  the  fatigue 
performance  of  structural  materials.  This  present  study  further  develops  one  phase  of 
that  work  on  applying  reliability  methodology  to  the  fatigue  performance  of  aircraft 
structures. 

At  best,  fatigue  behavior  is  a complex  one.  It  concerns  not  only  the  spectrum  of 
variables  like  loads,  geometry,  structural  response,  manufacture,  etc.,  all  affecting  the 
initiation  and  extension  of  fatigue  damage,  (i.e.,  the  fatigue  crack)  but  also  the 
coincident  material  behavior.  Of  particular  note  is  the  variability  in  the  response  of 
materials  to  this  complicated  total  environment.  Scatter  in  the  fatigue  performance  of 
materials,  exposed  to  a.  given  load  environment,  is  an  observation  well-documented  in 
the  annals  of  structural  mechanics  ever  since  fatigue  behavior  was  first  identified.  The 
early  unpredicted  appearance  of  fatigue  damage  as  a detectable  crack,  iol lowed  by  its 
subsequent  propagation  to  cause  failure  if  undetected,  has  a potentially  detrimontal 
effect  on  the  function  and  safety  of  operational  structures,  not  to  mention  the  serious 
economic  consequences.  However,  with  the  development  of  the  science  of  fracture 
mechanics  during  the  last  20  years,  a means  has  evolved  by  which  a structural 
material’s  fatigue  performance  can  be  characterized  for  the  conditions  ol  crack 
propagation  and  residual  strength. 

Obviously,  the  measured  time  to  fatigue  crack  initiation  in  a structural  system  is  a 
function  of  the  detection  capability  for  the  crack  us  well  as  the  complicated  cumulative 
response  to  its  environment.  The  crack  growth  rate  and  fracture  toughness  of  a material 
are  physical  parameters  identifying  the  other  important  proportios  necessary  to  assess 
the  durability  and  integrity  of  a structure  during  its  operational  lifetime.  These  throe 
fatigue-related  properties  exhibit  variability  and  depend  on  experimental 
determination.  Accordingly,  this  empiric  nature  of  the  fatigue  process  makes  it  a ready 
candidate  for  application  of  statistical  techniques  and  reliability  technology. 

As  a mutter  of  coincident  interest,  tho  strength  properties  of  materials  are  treated  with 
a well  tested  and  widely  accepted  reliability  methodology  in  defining  the  mechanical 
properties  of  strength.  A normal  distribution  is  assumed  for  a particular  property  like 
ultimate  tensile  strength.  This  mechanical  property  is  measured  and  distribution 
parameters  estimated  from  the  data,  allowing  a reliable  level  of  strength  to  bo 
established.  Thus,  sampling  techniques  and  the  selection  of  some  degree  of  assurance 
provide  a confidence  and  probability  level  with  which  the  spooific  alloy  and  condition 
will  demonstrate  the  desired  strength  property. 

However,  the  task  of  applying  the  mutorinl  to  a structural  configuration  introduces 
some  additional  difficultios.  Major  portions  of  the  structural  system  fulfill  a multiple 
role  of  providing  aerodynamic  properties  and  function  and  structural  capability  at  very 
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stringent  levels  of  structural  efficiency.  This  situation  can  introduce  some  hazards  in 
the  design  process.  Nevertheless,  ultimate  <itrength  goals  usually  present  only  minimal 
difficulties  in  matching  a prescribed  load  envelope  to  a structural  capability.  A 
well-developed  analysis  technology  backed  by  an  experimental  strength  verification 
procedure  of  static  testing  generally  assures  achieving  a specified  strength  goal. 
Structural  strength  deficiencies  uncovered  during  this  process  are  usually  resolvable  by 
effective  modifications,  often  simple,  to  achieve  the  specified  strength  levels  within  the 
production  fleet. 

On  the  other  hand,  attaining  durability  is  a moro  difficult  task.  It  has  a more  complex 
design  and  verification  process,  since  \ is  subject  to  unexpected  and  unidentified 
operational  exposure  than  that  presumed  during  design.  Furthermore,  the  analysis 
process  apparently  seems  neither  as  wefi  developed  nor  as  universally  understood  as  the 
comparable  static  strength  design  procedures.  Likewise,  remedial  steps  are  not  as 
simple  as  those  for  static  strength  because  of  the  cumulative  aspects  rather  than  the 
single  discrete  response  to  a load.  A simple  corrective  effort  may  introduce  other  fatigue 
sensitive  details  and  thus  be  only  temporary  in  nature.  Additionally,  fatigue  damage 
initiation  and  propagation  are  continuing  hazards  over  the  total  operational  life  of  the 
structural  system.  Hence,  the  capability  for  analysis  must  not  only  encompass  the 
Initial  paper-to-hardware  development  phase  but  also  need  to  be  capable  of  coping  with 
results  from  the  operational  phase.  However,  despite  the  very  local  nature  of  fatigue 
damage  and  its  difficulty  of  detection  in  structural  details,  some  degree  of  damage 
tolerance  built  into  the  surrounding  area  can  provide  the  opportunity  to  extend  the 
serviceability  of  the  structure  beyond  the  first  bell-ringer  incident  of  fatigue  crack 
initiation. 

Approaches  compensating  for  fatigue  life  variability  in  aircraft  structural  systems  have 
ranged  from  selecting  arbitrary  scatter  factors,  like  2 or  4,  to  using  computational 
procedures  based  on  theory.  As  in  the  ultimate  strength  problem,  the  primary  emphasis 
on  developing  a rationele  to  characterize  the  scatter  in  fatigue  behavior  has  been  to  fit 
a statistical  distribution  function  to  observed  data.  With  an  identified  distribution 
function  and  estimated  parameters,  the  likely  performance  and  the  associated  factor  of 
safety  to  assure  an  acceptable  level  of  fatigue  performance  are  determinable.  The 
measure  of  success  of  this  approach  is  the  accuracy  with  which  the  selected  distribution 
and  its  parameters  replicate  the  observed  scatter  in  fatigue  performance.  The  particular 
advantage  (i.e.,  gaining  structural  efficiency)  in  this  approach  is  the  reduction  in 
tolerance  gained  through  identifying  only  the  shape  and  scale  parameters  of  the 
distribution  function.  For  each  structural  detail,  the  determination  of  the  location 
parameter  is  always  necessary.  Added  assurance  is  gained  on  the  likely  performance  as 
more  data  are  acquired  to  estimate  the  location  parameter  (i.e.,  the  mean  or  other 
characteristic  life).  Then  the  probabilistically  determined  scatter  factor  applied  to 
characteristic  life  provides  a practical  way  to  attain  reliability  in  fatigue. 

There  are  several  problems  associated  with  this  approach  to  the  fatigue  damage 
initiation  prediction  previously  described.  B'irst,  identification  of  the  distribution  is  done 
empirically,  although  the  central  tendency  (the  mean,  median,  or  characteristic  life)  is 
reasonably  identified  by  a few  test  values.  However,  for  a largo  group  of  identical 
structures  in  a specific  operating  environment,  the  operational  condition  of  the  group  is 


more  closely  identified  by  the  time-to-firrit  failure  (or  subsequent  early  failures)  of  a 
detail  rather  than  the  average  behavior.  Consequently,  to  determine  the  adequacy  of 
the  fit  of  the  distribution,  the  quantity  of  test  data  must  be  sufficient  to  verify  the 
accuracy  of  the  toe  of  the  distribution.  Accordingly,  hundreds  or  even  thousands  of 
detail  fatigue  test  results  are  necessary  to  recognize  and  estimate  the  probability  of  the 
early  failure  in  the  group  of  parts.  Collecting  such  extensive  data  may  be  practically 
and  economically  impossible  for  a particular  structure.  Another  facet  of  the  problem 
concerns  the  very  nature  of  a statistical  behavior;  the  actual  life  of  a part  is  not 
predictable  on  a part-to-part  basis  but  instead  only  by  the  probable  performance  as 
determined  by  a general  distribution.  Thus,  the  hazards  of  population  identification 
from  samples  must  be  resolved  directly  after  incidents  occur.  Alternatively,  the  parts 
may  be  replaced  on  the  basis  of  an  acceptable  probability  of  failure  from  a safe  life 
reliability  analysis  scheme.  Of  course  this  makes  some  assumptions  about  the  expected 
usage  of  a part.  However,  a fail-safe  or  damage-tolerant  structural  design  process 
provides  a practical  opportunity  to  treat  fatigue  crack  initiation  on  a part-to-part  basis 
if  unanticipated  damage  does  initiate  prior  to  predicted  goals. 

Some  of  the  early  work  sponsored  by  AFML  explored  the  variability  of  the  initiation  of 
fatigue  damage  in  terms  of  statistical  and  related  reliability  concepts.  This  led  to  an 
examination  of  the  time-to-first  failure  rather  than  a factored  "average”  time  to  assess 
fatigue  performance.  In  reference  1,  characteristics  of  the  log-normal,  Weibull,  and 
gamma  distributions  were  studied  relative  to  the  first  failure  for  a range  of  sample  (i.e., 
fleet)  sizes.  A further  development  on  the  application  of  this  first  failure  concept  to  a 
fleet  was  reported  in  reference  2.  This  latter  investigation  developed  the  concept 
involving  specific  structural  details  in  relation  to  the  task  of  design,  manufacture,  and 
subsequent  operational  exposure.  This  study  considered  the  associated  practices  of 
fatigue  life  estimation,  supported  by  developmental  and  verification  testing  of  the 
aircraft  structural  system.  Primary  emphasis  focused  on  the  application  of  the 
log-normal  and  Weibull  distributions,  describing  fatigue  damage  initiation  variability. 
Scale  parameters  were  determined  by  examining  tost  data  from  many  groups,  mostly 
small  ones,  of  various  types  of  specimens,  A standard  deviation.of  0,14  was  derived  for 
the  log-normal  distribution  shape  parameter,  while  the  value  4.0  was  appropriate . for 
the  Weibull  distribution.  Using  the  Weibull  distribution  function,  an  application  of  the 
time-to-first  failure  concept  was  made  to  a tanker  transport  type  of  aircraft  subdivided 
by  usage  into  a number  of  different  fleets.  Eight  different  details,  which  had  laboratory 
identified  fatigue  performance  during  the  design  stage,  were  investigated  for  their 
subsequent  behavior  in  the  operational  fleets.  A promising  correlation  was  found 
between  the  predicted  fatigue  performance  and  service  results  available  at  that  time. 

In  reference  3,  the  problem  of  identifying  a representative  distribution  to  characterize 
fatigue  scatter  was  explored  for  aluminum  alloy  materials  as  well  as  titanium  and  steel. 
A technique  for  normalizing  and  utilizing  the  data  from  many  groups  of  tests  was 
applied.  This  process  made  available  thousands  of  test  points  and  provided  an  insight 
into  the  extreme  behavior  of  the  normalized  test  data.  The  Weibull  distribution 
appeared  to  fit  the  tendency  of  the  data  to  flare  out  at  very  low  probabilities.  Other 
subsequent  studies  were  made  on  the  application  of  this  reliability  technology  scheme  to 
another  type  of  cargo  airplane  and  a fighter  aircraft,  references  4 and  5,  respectively. 
Both  of  these  studies  indicated  that  the  early  failure  concept  provided  a corroboration 


with  fleet  performance.  However,  the  Weibull  distribution  provided  a more  conservative 
estimate  of  fatigue  damage  initiation  than  the  log-normal  distribution.  Furthermore, 
the  fatigue  results  demonstrated  by  laboratory-programmed  loading  tests  indicated 
scale  factors  less  than  the  standardized  values  of  reference  2.  This  is  not  a 
contradiction,  since  it  is  the  result  of  statistical  assessment  of  observed  data  and  the 
related  problem  of  sampling. 

In  the  reference  6 study,  two  other  concepts  were  explored  for  application  of  the 
reliability  analysis  scheme  utilizing  the  first  failure  approach.  One  proposed 
determination  of  the  scatter  factor  as  the. ratio  of  two  random  quantities,  namely  the 
test  life  and  service  life  from  a Weibull  population.  This  approach  is  different  from  the 
previously  proposed  scheme  which  samples  a population  to  estimate  its  location 
parameter  and  calculates  the  likely  occurrence  of  the  first  or  early  failure.  The  smaller 
scatter  factor  of  the  former  estimation  scheme  is  illustrative  of  the  different  approach. 
A second  variation  was  the  definition  of  a joint  scatter  factor  to  combine  the  effects  of 
fatigue  and  usage.  In  the  reference  2 approach,  the  characteristic  or  mean  life  is  treated 
as  an  input,  principally  because  of  the  very  nature  of  cumulative  fatigue  damage 
response  of  materials.  It  is  of  pertinent  interest  to  note  that  both  references  4 and  5 had 
illustrations  in  which  prediction  did  not  correlate  with  service  results.  However,  this 
observation  is  not  unique  to  those  two  investigations.  Hence,  the  introduction  of  a joint 
scatter  factor  for  fatigue  and  use  is  not  included  in  the  basic  reliability  analysis  scheme. 
Certainly,  additional  work  seems  necessary  to  evaluate  and  judge  the  part  played  by 
measured  usage  in  terms  of  cumulative  acceleration  counts  and  the  responding 
cumulative  damage, 

Considering  further  the  choice  of  a representative  distribution  to  simulate  the 
variability  in  time-to-fatigue  crack  initiation  under  a given  loading  environment,  an 
experimental  study,  summarized  in  references  6,  7,  and  8,  generated  order-statistic  test 
data  from  multideteil  specimens.  A graphical  study  of  all  the  data  (refs.  7 and  8)  did  not 
indicate  an  obviously  superior  fit  of  the  experimental  external  data  by  either  the 
Weibull  or  log-normal  form  of  distribution.  However,  the  Weibull  distribution  generally 
predicted  an  earlier,  more  conservative  estimate  of  the  first  failure  than  did  the 
log-normal.  This  was  particularly  apparent  as  the  level  of  reliability  varied  from  0.60  to 
0.96  for  the  early  failures. 

As  noted  in  a preceding  paragraph,  the  log-normal  and  two-parameter  Weibull 
statistical  distributions  were  compared  with  empirical  distributions  obtained  from  a 
considerable  amount  of  normalized  fatigue  test  data  (ref.  3).  These  empirical 
distributions  were  observed  to  have  a tendency  toward  symmetry  between  the  upper  and 
lower  halves  when  plotted  on  normal  probability  paper,  and  also  to  have  a flexure  mode. 
This  latter  behavior  disagrees  with  the  charact.eristicB  of  the  log-normal  distribution, 
whereas  the  former  symmetry  condition  is  not  found  in  results  generated  from  a 
Weibull  model.  However  by  transforming  the  basic  two-parameter  Weibull  distribution, 
a new  "symmetric-Weibull"  distribution  was  obtained,  which  had  the  baoic  form  of  the 
empirically  obtained  fatigue  test  distribution. 

Subsequent  to  this  investigation,  it  was  determined  that  the  modified  or 
symmetric-Weibull  distribution  belonged  to  a new  class  of  statistical  distributions 
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(ref.  9)  that  contained  certain  reciprocal  properties  (ref.  10).  It  was  clear,  therefore,  that 
other  distributions  within  this  class  should  be  examined  to  determine  if  a more  suitable 
fatigue  life  distribution  model  was  available.  In  this  context  the  term  suitable  is  not 
necessarily  limited  to  an  evaluation  based  on  data  correlation  but  must  include  aspects 
of  ease  of  application,  philosophic  rationale,  etc.  For  example,  the  fatigue  distribution 
model  that  would  show  the  most  correlation  with  the  data  is  one  formed  empirically 
from  the  data.  This  would  contain  many  more  parameters  (ref.  11)  and  would  thus 
beOome  extremely  difficult  if  not  impossible  to  apply  to  typical  problems  such  as 
parameter  estimation  from  censored  observations,  determination  of  confidence  bounds, 
and  so  forth. 

It  should  be  realized  that  an  airplane  designer  is  frequently  guided  by  the  past 
performance,  either  operational  or  test,  of  similar  components  to  the  one  presently  being 
designed  or  modiiled.  Operational  fatigue  data  are,  without  exception,  of  the  extreme 
event  t^rpe,  and  as  such  information  needs  to  be  assessed  it  is  necessary  to  have  an 
estimation  procedure,  preferably  giving  unbiased  values,  for  the  general  case  of 
censored  data.  Consequently  it  is  advisable,  for  reasons  of  practicality,  to  keep  the 
numbers  of  parameters  or  variables  to  a minimum. 

The  correct  interpretation  of  operational  data  is  further  confounded  by  the  fact  that  ail 
of  the  limited  information  may  not  bn  completely  identical,  but  come  i>om  nominally 
different  subsets;  e.g.,  small  differences  in  failure  locations,  structural  operation,  etc.  It 
is,  therefore,  most  desirable  to  have  the  capability  of;  (1)  obtaining  bounded  estimates  of 
the  distribution  parameters  iVom  the  censored  data  and  (2.)  testing  the  hypothesis  that 
different  censored  data  subsets  actually  do  or  don't  come  from  the  same  patent 
population. 

While  fatigue  crack  initiation  is  a key  milestone  in  the  usable  life  of  a structural 
£iystem,  the  presence  of  fail-safe  or  damage-tolerant  capabilities  in  structures  introduces 
concepts  of  crack  growth,  inspection  control  for  detection  of  initiated  cracks  before  the 
(iL'lerioration  of  structural  safety,  and  the  upper  limit  of  crack  growth  as  dictated  by  the 
residual  strength  capability  of  the  fatigue-damaged  structure.  These  elements  of  fatigue 
performance  of  a structural  system  are  all  a part  of  the  basic  reliability  plan  of  this 
current  investigation.  This  multidiscipiine  structural  reliability  method,  initially 
described  in  reference  12  and  subsequently  revised  and  expanded  in  references  13 
and  14,  recognizes  that  the  time-to-crack  initiation  and  its  growth  id  a random  event. 
However,  current  Air  Force  philosophy  (ref.  16)  retlects  an  approach  based  on  the 
preexistence  of  a flaw  at  a critical  detail  in  the  structure  when  it  enters  operational 
service.  Hence,  some  accounting  for  the  concept  of  an  undetected  flaw  at  time  zero  is 
also  a desirable  element  of  this  reliability  analysis  scheme. 

This  report  discusses  the  findings  from  an  investigation  of  the  topics  mentioned  in  the 
preceding  paragraphs.  The  new  class  of  distributions  is  discussed  in  section  2.0,  and  four 
distributions  from  this  class  are  highlighted.  Section  3.0  discusses  the  development  of 
parameter  estimation  procedures,  both  censored  and  uncensored  data,  for  one  of  the 
described  distributions.  An  efficient  method,  which  is  believed  to  be  a major  advance,  of 
obtaining  unbiased  estimates  and  confidence  bounds  for  censored  data  is  also  described 
in  this  section.  Additionally,  a means  of  testing  the  hypothesis  that  two  independent 
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i;  data  samples  belong  to  the  same  parent  population  is  presented.  Finally  these  developed 

procedures  have  been  applied  for  illustrative  purposes  to  some  real  samples  of  censored 
fatigue  data. 

; The  discussion  given  in  section  4,0  of  this  report  covers  the  rationale  for  development  of 

; the  preflaw  model;  section  5.0  describes  the  model  and  its  application  to  the  existing 

reliability  analysis  system  given  in  references  12  and  13.  Finally,  section  6.0  covers  the 
results  of  some  parametric  studies  made  using  the  expanded  reliability  system. 

I'- 

Two  Other  sections  complete  the  document.  Section  7.0  discusses  structural  data  input 
considerations,  and  section  8.0  gives  the  conclusions  reached  from  this  investigation  and 
lists  some  recommendations. 
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2.0  THE  ADVANCED  DISTRIBUTIONS 


2.1  A NEW  CLASS  OF  ADVANCED  DISTRIBUTIONS  FOR  FATIGUE  DATA 

A considerable  amount  of  fatigue  test  data  was  collected  and  normalized  in  an  earlier 
investigation  (ref.  3).  These  normalized  results  were  then  used  to  define  empirical 
distributions  which,  when  plotted  on  normal  probability  paper,  demonstrated  both  a 
flexurs  at  the  tails  and  a symmetry  between  upper  and  lower  halves. 

A transformation  of  the  two-parameter  Weibull  into  a symmetrical-Weibull  distribution 
was  developed  and  was  found  to  have  thd  characteristics  of  the  einpirical  distribution. 
Subsequently,  it  was  established  that  the  symmetric-Weibull  distribution  belonged  to  a 
new  advanced  class  of  distributions  that  had  certain  reciprocal  properties  (refs.  9 
and  10). 

The  following  initial  study  considers  four  distribution  ihnctions  from  within  this  new 
class  of  distributions.  The  first  is  the  log-normal  distribution  that  can  be  used  as  the 
reference  baseline,  and  another  for  obvious  reasons  is  the  symmetric-Weibull.  Two  other 
distributions  are  those  from  reference  10  and  are  labelled  as  the  hyperbolic  sine  and  the 
inverse  hyperbolic  sine  distributions. 

These  distributions,  described  in  the  following  paragraphs,  are  based  on  three 
parameters.  In  their  application  to  modeling  fatigue  performance,  it  is  assumed  that  a 
flexure  parameter  will  represent  material  behavior;  u scale  parameter  will  reflect  such 
' variables  as  configuration,  size,  type  of  loading,  etc.,  which  are  characteristics  of  the 
specimen  and  the  fatigue  test;  and  finally  a location  parameter  which  will  depend 
primarily  on  the  magnitude  of  the  loads  to  which  the  test  specimen  is  subjected. 

2.1.1  THE  f-NORMAL  DISTRIBUTION 

Let  5 be  a class  of  real  valued  functions  so  that  if  then  it  sat  , sues  the  functional 
equation 


t(t)  = -t(l/t)  fort>0  (2.1.M) 

If  T is  a nonnegatlvo  random  variable  so  that  there  exists  constants  a,/3  > 0 and  (tB, 
for  which 


iuj)  ~ Normal  (0,1) 


then  by  adopting  a nom^'nclature  analogous  to  that  of  the  log-normal  distribution,  it  can 
be  said  that  T is  ^-normal. 
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Thus  if  T is  ^-normal  with  parameters  a,li>0  and  has  a distribution  F then 

\ 


F(t)  = 4>[~5(t/^)]  fort  > 0 


■;  I 


where 


4>(x)  = f 
‘s/2ir  J 
00 


(2.1.1-2) 


To  facilitate  comparison  with  currently  available  plotted  data  that  are  so  frequently  in 
terms  of  the  logarithm  of  the  number  of  cycles,  consider  the  log>life  variate  X = 8nT  and 
take  its  distribution  G to  be  of  the  form 


G(x) 


<I»[-^co(x-M)J 


-«>  < X < “O 


(2.1. 1-3) 


where  = fai/3;  o)(x)  =•  ^(o*)  for  all  x. 

Now  consider  a class  of  functions  SI  defined  as  follows: 

ioefi  iff  w'  > 0,  oj'  > oi"  and  o>  is  odd,  checking  that  weO  iff  fceE 

whencc(x)  * -«i  < x < «» 

Examining  now  the  closure  properties  of  SI.  Assuming  that  if 


(2.1.1-4) 


and  moreover  if 


cj, , cOj  then  w,  + cOj  eSl 


ioeSl  and  a > 0, 0 < b < 1 


then  both 


co,,Wj  eS2  where  CO,  (x)  = aco(x)  and  Wj(x)  = co(bx)  for-“»  < x < «> 


As  a consequence,  for  any  given  wefl,  a three-dimensional  parametric  family  of 
distributions  of  the  log-life  can  be  generated,  namely. 


(2.1.1-5) 


G(x:  a, 7^)  *=  -«>  < x < <» 

where  «>0,  y>0,  ~ou<^<oo. 

In  thia  notation,  n is  the  location  parameter,  y is  the  scale  parameter,  and  a is  the 
flexure  parameter.  The  flexure  and  scale  both  together  control  the  shape.  (Note  that  if 
it  is  required  that  <a' >m'\  a restriction  y > 1 must  apply  to  ensure  the  closure 
properties.) 

Examples  of  such  functions  a,  which  correspond  to  known  parametric  families,  for 
-»  < X < w,  are: 

(0,(x)  s sinh  X 

Wj(x)  = <o,"*(x)  = en(x+'/x^T) 

which  corresponds  to  the  distributions  of  Birmbaum-Ssunders,  and 

WjCx)  “ d>'[F(x)] 

where 

F(x)  = + (sgn  X)  exp  {-(Cn2)e‘'^'| 

which  corresponds  to  the  symmetric-Woibull  distribution. 

Finally  for  comparison,  a well-known  and  commonly  used  fatigue  life  distribution  is 
taken  as 

= X 

which  corresponds  with  the  log-normal  distribution  as  it  was  shown  earlier  that 
X*finT. 

2.1^  DISCUSSION  OF  EXPLORATORY  APPLICATION  AND  RESULTS 

The  typical  behavior  of  this  class  of  distributions  is  illustrated  in  Hgures  1 through  4. 
These  show  cumulative  frequency  plots  of  the  four  statistical  models,  discussed 
previously,  for  several  arbitrary  values  of  the  flexure  parameter  a.  It  is  clear  from  these 
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illuBtrations  that  the  log-normal  model  (fig.  1)  plots  linearly  and  thus  does  not  have  the 
flexure  which  has  been  observed  to  be  characteristic  of  fatigue  data.  The  other  three 
selected  models  are  all  seen  to  flex  with  the  hyperbolic  sine  diatrlbution,  (fig.  2),  flexing 
in  a direction  opposite  to  either  the  inverse  hyperbolic  sine  or  the  symmetric-Weibull 
distributions,  figures  S and  4,  respectively.  These  latter  three  figures  also  vividly 
demonstrate  the  importance  of  the  flexure  parameter  a,  because  it  is  obvious  how  the 
frequency  plot  can  go  from  almost  linear  to  highly  flexed  depending  upon  the  magnitude 
of  a.  Finally  the  scale  parameter  y is  seen  to  govern  the  slope  of  the  plots  and  is  thus 
analogous  to  the  standard  deviation  parameter  of  a normal  distribution. 

A large  amount  of  aluminum  fatigue  data  was  normalized  and  plotted  on  normal 
probability  in  reference  3.  These  data  have  been  reproduced  in  this  report  in  figure  6 
using  the  a-ordinate,  y-abscissa  format  of  figures  1 through  4.  The  amount  of  flexure 
exhibited  by  the  data  is  immediately  obvious,  and  the  almost  "mirror  image"  nature  of 
the  lower  and  upper  halves  of  the  curve  is  also  clear.  The  best  log-normal  and 
symmetric-Weibull  distribution  values  are  superimposed  on  the  data-curve.  This  chart 
demonstrates  that  the  log-normal  plots  linearly  as  expected  and  the  symmetric-Weibull 
does  follow  the  data  trend  and  flex,  although  not  r gh  to  exactly  match  the  data  at 
the  extreme  tails. 

Figure  6 shows  the  comparison  of  the  Birmbaum-Saunders  inverse  hyperbolic  sine 
distribution  with  the  aluminum  data  curve.  A value  of  a » 0.74  was  used  to  compute  the 
curve  shown  in  this  figure.  It  is  clear  from  figure  6 that  this  distribution  model  can 
virtually  duplicate  the  empirical  distribution,  with  only  a small  minmatch  in  the 
extreme  upper  tail  at  approximately  the  upper  two-  or  three-tenths  of  1% . 


It  had  been  presumed  at  the  beginning  of  this  investigation  that  the  flexure 
parameter  a would  bo  a material  parameter  and  that  the  scale  parameter  y would 
reflect  stochastic  differences  resulting  from  variations  in  specimen  geometry,  testing 
methods,  etc.  Thus,  it  can  be  anticipated  that  the  flexure  parameter'  found  to  produce  a 
good  fit  with  all  the  aluminum  data  will  bo  found  to  bo  close  to  the  individual  values 
obtained  after  analysis  of  the  eeparate  data  subsamples.  Furthermore,  in  order  to 
facilitate  the  development  of  the  subsequent  estimation  theory  and  its  applicutioit,  it 
was  decided  at  this  time  to  concentrate  on  the  inverse  hyperbolic  sine  distribution. 
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3.0  DISTRIBUTION  PARAMETER  ESTIMATION 

3.1  THE  ESTIMATION  OF  PARAMETERS 

Suppose  that  there  is  a complete  sample  of  failure  data  (xi,.,.,xn)  from  the 
distribution 


G(x)  = 


oo  < \ < OO 


where  m is  known  but  the  parameters  a,  y,  ii  are  unknown.  The  density  of  each 
observation  is 


and  thus 


fing(x)  “ -ftn(a7)+Cnw'^~^^ 

The  log-likelihood  then  is 


L **]C*2ng(Xj') 

which,  except  for  a constant  independent  of  the  parameters,  is 


n 

L =2: 
1 


■fin(ory)-~w* 
2a* 


In  order  to  maximize  the  likelihood  consider 


at 

3p 


3L 


1 VXi-M\ 

^ It  / 

/Xi-P\ 

^\y  ) 

\ / 

lk»V 
0a  ~ 4^ 


a 


(3.1-1) 


(3.1-2) 


(3.1-3) 


(3.1-4) 


Now  the  joint  solution  of  dhldii'=‘0,  dLldy^O,  dh/da  -^0  will  yield  the  maximum 
likelihood  (ML)  estimators  cf,  y,  and  /x.  Thus,  the  simultaneous  solution  to  the  following 
equations  in  the  variables  (o,  y,  /x)  is  sought; 

)=E-^  (3.1-5) 

(3.1-6) 


As  a notational  convenience,  define  as  Implicit  iVinctions  of  /x  and  y the  values 


y.  s -y  i = 1 , . . n 


and  introduce,  for  any  function  f,  the  averaging  operator 


<iXyi)>  =”Ef(yi) 


r 


Define  the  functions  of  (/x,  >i: 


" = SI  ■ 


H*  = (yiR{yi))-l 

<coMyi)> 

w .ere  the  ratio  and  product  functions  are  defined  by; 

R =>  w"/w',  P = CjUw' 


and  both  R and  P are  odd. 


3.1.1  MAXIMUM  LIKELIHOOD  ESTIMATE  (MLE)  ITERATION  METHOD 

The  iteration  method  for  the  situation  when  three  parameters  are  unknown  is  as  follows.  For 

0,1,2,... 

Step  1:  Given  /nj,  solve  for  y in  the  equation 

H*  (/ij,7)  = 0 

Call  the  solution  yj 

Step  2:  Given  yj  bolve  for  /t  in  the  equation 

H(M,7j)  “ 0 


Call  the  solution  pj+i 
Set 


then  set 


and  ju  ■=  Hm  jui,  ^ = liin  71 


-1- 


where 


^f(x^  is  the  averaging  operator  = 


To  apply  the  Newton-Raphson  technique  above  to  solve  the  equation  Rx)  »0  for  the 
single  variable  set 

f(xi) 

Xj+1  “ Xi  ' — — j “ 0,  1 , 2, . . . 


then  lim  xj  xq  will  satisfy  f)(xo)  "*0. 

j-b<» 
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Thus,  it  is  necessary  to  know  both  partial  derivatives 


07  ’0/11 


so  as  to  be  able  to  determine  their  behavior. 
Examining  first  dH/dfi  ; it  is  checked  that 
9 /«/..  A _ 1 


•^(P(yi))  = -;J-<P'(yi)> 
(R(yi)>  - (R'(yi)> 


(3.1.M) 


9H  - 


<p(yi)> 

(w^Cvi))  * 


■+(R'(yi)> 


Similarly  it  is  found  that 
0H' 

^57 


. ^ H>i)> 


^ <"  O'!*)  <w>(y|)> 


(3.1.1.2) 


3.1.3  DISCLSSION  OP  EXPLORATORY  APPLICATION  AND  RESULTS 

The  initial  task  in  testing  the  MLE  equations  and  the  Newton>Raphson  iterative 
techniquee  (sec.  3,1.1)  was  the  selection  of  some  suitable  data.  A large  number  of  such 
samples  are  listed  in  appendix  III  of  reference  2.  In  order  to  facilitate  the  initial 
estimates,  it  was  decided  to  select  data  of  sample  size  s>  10  specimens.  The  first  such 
acceptable  uncensored  test  sample  comprising  10  specimens  was  selected  from 
reference  2,  and  the  described  MLE  methodology  was  applied.  Unfortunately  at  this 
initial  attempt  at  estimating  the  distribution  parameters,  it  was  found  that  the 
estimators  .did  not  exhibit  convergence  and  so  no  solution  was  obtained. 
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An  understanding  of  the  nature  of  this  problem  was  most  important,  and  so  it  was 
decided  to  generate  several  values  of  the  partial  derivatives  dhld/x  - 0 and  dUdy  ==  0, 
(see  sec.  3.1).  The  plotted  results  shown  in  figure?  demonstrate  that  both  quantities 
appear  to  be  approaching  their  respective  asymptotic  values  without  ever  intersecting; 
i.e.,  no  common  value  as  would  be  necessary  for  an  MLE  solution,  Thus  it  seems  that  for 
some  quite  typical  fatigue  data,  as  evidenced  by  the  test  results  given  in  figure  7, 
problems  would  be  encountered  using  the  present  form  of  the  MLE  equations. 


To  obviate  the  difficulty  encountered  in  this  numerical  analysis,  equivalent  forms  of  the 
MLE  equations  were  tried  but  without  any  more  success.  Attempts  at  redefining  or 
reparameterizing  the  equations  also  proved  fruitless.  Subsequent  additional  studies 
then  showed  that  in  certain  cases  the  simultaneous  MLE*s  of  all  three  parameters  do 
not  exist. 

3.2  NONEXISTENCE  OF  THREE  MLE’S  IN  CERTAIN  CASES 


In  section  3.1.2,  mention  was  made  of  the  lack  of  an  MLE  solution  for  a fairly  typical 
sample  of  fatigue  data.  The  data  values  and  the  estimator  results  wore  shown  plotted  in 
figure  7.  Consider  how  the  behavior  of  the  likelihood  function  whon  is  S-uhaped  ns 
shown  in  the  schematic  below; 


Now  0)  is  S'shaped  if 


1 > > co'(x)  > 0 w'(0)  “ 1 


co'(x)  and  — ^ — both  decrcuso  as  1x1  increases 


(3.2-2) 


Cotiaider  the  joint  log-likelihood  equation  involving  only  the  parameters  /a,  y having 
previously  eliminated  a.Then  write 


L(m,7)  = -8n(onr) 

2na* 


where  a is  a function  of  /u.,  y given  by 


After  substitution  and  simplification  it  is  found  that 


(3.2-8) 


Now  set 


Vj  = Xj-/i 

i = 1, ....  n, 
n 

a = nvix  |yi| 
i=l 


Thus  from  equation  (3.2-1)  for  i»l,  . . n 

[ycj(yi/7)]^  > yj*[w'(yi/7)]^  > [w'(a/7)]  ^ 
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Upon  substitution  it  is  found  that 


Uti.y)  + 


nZIj  Vi^f  + 2 - £nw'(a/7) 


On  the  other  hand  from  the  bound  in  equation  (3.2-1) 


[7w(Vi/7)V  = < Vi*  ‘ J . • • •.  n 


end  thus 


Mm, 7)  -5+^Ej2nw'(yj/7) 


From  equation  (3.2-2)  it  is  known  that 

<^'(y{h)  > w*(a/7>  for  i “ 1 , . . n 


and  thus  from  equation  (3.2-4)  it  follows 


L(M,7)  +2*^n 


nZ)j  yi*['^2  ^ *2nw'(a/7) 


As  a result  it  may  be  remarked  that  the  log-likelihood  is  bounded  by: 

8nw'(a/7)  < Um,7)  (xj-M)*j  < -finw'(a/7) 


where 

n 

a = max  |x[-/i| 
i=l 


(3.2-4) 


It  follows,  of  course,  that  in  the  limit  as  the  log-likeUhood  approaches  a limit  for 
which  the  maximum  value  of  ii  is  at  3C.  Moreover  those  (o’s  chosen  for  study  show  that 
this  limit  is  reached  very  quickly.  In  addition  computation  also  indicates  that  with  ui 
defined  as  shown  earlier,  that  with  fixed,  L(fi,y)  is  an  increasing  function  of  y > 0. 

Thus  the  likelihood  is  maximum  at  <=  x,  and  x is  an  unbiased  estimate  of  fi.  But 

the  maximum  likelihood  of  y does  not  exist. 

The  problem  of  characterizing  the  class  of  a)’s  for  which  U/t.y)  is  an  increasing  fuhction 
of  y seems  to  be  difficult.  It  must  be  that  as  y*»  the  second  term  of  equation  (3.2-3) 
increases  while  the  first  term,  which  decreases,  is  completely  dominated. 

Faced  with  the  fact  that  for  some  a>'s,  which  are  of  practical  interest,  the  MLE’s  of  all  three 
parameters  may  not  exist,  it  is  now  necessary  to  seek  some  alternate  method  of  estimating 
the  flexure  parameter.  This  follows  in  the  next  section. 

3.3  A SIMPLIFIED  ESTIMATOR  FOR  UNCENSORED  DATA 

The  discovery  that,  for  some  fairly  typical  samples  of  data,  MLE’s  of  all  three 
parameters  do  not  exist  meant  a revision  to  the  planned  approach.  Now  it  was  believed 
that  if  one  of  the  throe  parameters  could  be  assumed  to  have  a known  fixed  value,  then 
MLE’s  could  be  made  of  the  remaining  two  parameters.  The  most  likely  candidate  for 
having  least  variation  was  considered  to  be  the  flexure  parameter  a,  and  so  some  means 
were  necessary  to  determine  its  probable  value. 


In  order  to  facilitate  this  development  stage  of  the  investigation,  it  v/as  decided  to 
concentrate  on  the  fairly  large  uncensored  samples  of  data  from  reference  2 in  the 
determination  of  the  distribution  parameters.  Keeping  the  problems  of  the  MLE  in 
mind,  some  consideration  was  given  to  the  development  of  moment  estimators  as  a 
simpler  means  of  obtaining  parameter  estimates  from  fatigue  test  results. 
Unfortunately,  a general  method  suitable  for  all  distributions  within  the  class  of 
distributions  studied  was  not  determined.  Even  conventional  estimators  such  as  the 
method  of  least  squares  did  not  result  in  computable  estimates  of  all  the  location,  scale, 
and  flexure  parameters  simultaneously. 

However,  a fairly  simple  estimator,  ^hich  has  been  labelled  the  "STAR"  estimator,  was 
developed  by  combining  the  momeni.  and  MLE's.  The  application  of  this  estimator  was 
limited  to  the  Birmbaum-Saunders  inverse  hyperbolic  sine  distribution  and  has  shown 
itself  capable  of  consistently  solving  for  all  three  parameters  from  fairly  large  samples 
of  uncensored  fatigue  data.  For  example,  considering  the  distribution  given  as: 

o)^(x)  = sinh'*  (x)  (see  sec.  2.0) 
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and  the  solution  of  the  two  equations  in  a,  y;  namely, 


Now  solve  for  the  solution  y'*'  in  the  equation  in  y 


by  some  iteration  method.  Then 

It  is  known  that  fi*  is  unbiased  end  consistent,  and  so  o'**  and  y*  are  consistent 
estimators  and  have  been  found  to  converge  quite  readily. 

SJ.l  DISCUSSION  OF  EXPLORATORY  APPLICATION  AND  RESULTS 

Samples  of  fatigue  data  have  been  selected  from  appendix  III  of  reference  2 and 
analyzed  using  the  STAR  estimator.  Only  uncensored  samples  of  size  n » 6 specimens 
per  sample  were  selected.  Initially  it  had  been  the  intent  not  to  consider  any  samples 
comprising  leas  than  10  observations,  but  this  did  not  provide  a sufficiently  large  ^ta 
base,  and  so  additional  smaller  sized  samples  had  to  be  considered. 

The  results  of  the  data  analyses  have  been  plotted  in  the  form  of  cumulative  frequency 
versus  flexure  p.aramoter  graphs  (figs.  8 and  9).  It  must  be  noted  that  emphasis  has 
been  given  to  the  flexure  parameter  for  reasons  mentioned  previously,  although  all 
three  parameters  were  in  fact  computed  by  the  procedure  described  in  section  3.3. 

The  iS4  samples  of  fatigue  test  data,  which  were  collected  and  analyzed,  came  from  a 
multitude  of  original  sources  and  contained  several  variables  (table  1).  These  data  were 
consequently  grouped  in  several  different  categories  of  interest,  and  their  flexure 
parameter  characteristics  were  compared.  These  comparisons  are  illustrated  in  figures  8 
and  9.  It  is  readily  apparent  that  there  are  no  flindamental  differences  in  the  results 
presented.  All  show  that  the  large  mi^oi^ity  of  the  population  falls  within  a narrow 
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range  of  values  of  the  flexure  parameter.  This  small  variation  over  the  central  portion 
and  the  somewhat  larger  variations  at  the  extremities  are  most  probably  sampling 
error,  although  this  has  not  yet  been  established.  Furthermore,  it  is  noted  that  the 
weighted  mean  values  of  the  flexure  parameters  of  the  various  subgroups  are  all  very 
similar,  thereby  providing  additional  support  to  the  hypothesis  of  an  acceptably 
constant  value  for  this  parameter.  It  must  be  emphasized  that  the  data  selected  for  this 
study  covered  a considerable  range  in  the  extent  of  the  sample  scatter.  For  example, 
using  the  parameter  of  standard  deviation  of  logarithms  (base  10)  as  a measure  of 
variability,  the  largest  standard  deviation  within  the  analyzed  data  samples  had  a 
value  of  8 - 1.304.  However,  in  spite  of  gross  scatter  such  as  this,  the  flexure  parameter 
was  not  too  influenced,  thereby  demonstrating  its  relative  insensitivity  to  variability  in 
test  results. 

Finally  based  on  the  results  of  the  164  aluminum  alloy  test  samples  that  have  been 
analyzed,  the  weighted  average  value  of  the  flexure  parameter  is  a*  '=  0.64. 
Furthermore  having  determined  that  the  flexure  parameter  remains  almost  constant, 
the  three-parameter  distribution  model  being  investigated  can  be  treated  justifiably  as 
a model  with  one  known  constant  and  two  unknown  parameters  which  require 
estimating,  The  estimation  task  is  then  substantially  simplified  and,  therefore, 
facilitates  the  application  of  MLE  procedures  thus  providing  the  necessary  capability  of 
evaluating  smaller  sample  sizes  and,  more  importantly,  censored  data. 


3.4  ML  ESTIMATION  FOR  TWO  PARAMETERS  USING  CENSORED  DATA 

Assume  a data  sample  comprising  both  real  life  length  observations  X}, . . .,  x^  and  censored 
life  lengths  xi^^.  j , , . , , Xn.  In  these  circumstances  the  value  xi  is  taken  to  indicate  that  the  ith 
component  was  alive  at  time  xi  for  i » k+ 1 n. 

Assuming  that  the  life  length  observations  are  fbom  the  distribution 

G(x)  = d>  -‘»<X<oo 

where  u Is  a known  odd  function,  but  the  location  and  scale  parameters,  fi  and  y, 
respectively,  are  unknown.*  Then  the  density  of  each  observation  is 


*The  flexural  parameter  is  taken  an  being  constant  and  known  and  omitted  in  this  development 
by  assuming  a value  of  unity. 
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-and  flo.,,  ' 

fing(x)= ^Hny  + inw.'(^) 

The  probability  of  an  "alive-time**  of  value  x,  is 

P(X  > X)  = 1 - G(x)  = ♦ ^(^)] 

with  the  last  equality  following  since  u is  odd. 

It  follows  from  this  thht  the  log-likelihood  for  the  set. of  complete  data  (xi, . . Xn)  is 
Deflning  the  cumulative  hazard  as 

Q(x)  = -fin4»(-x) 


gives 


Upon  taking  partial  derivatives,  it  is  found  that 

_k 


9L  _ 1 1 , 


Cl!L)jTL)  1 "^\yJ\ 


i-k+l  L \ 7 /J  \ y h 


y 'y 


S-t|4-(=K.¥)fS)-^(¥)| 


■ 2L, +(?)]-<¥)(¥) 
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By  utilizing  the  product  and  ratio  functions,  P and  R,  respectively,  originally  defined 
in  section  3,1  as 


P “ ww'  and  R = co'Vco' 


and  setting 


h(x)»q{w(x)]w'(x) 


The  two  equations  dhibix  "•  0 and  dUdy  » 0 can  be  rewritten  as 


(3.4-1) 


and 


(3.4-2) 


Thus,  to  obtain  MLE*s  of  the  location  and  scale  parameters,  it  is  necessary  that 
equations  (3.4-1)  and  (3.4-2)  be  solved  simultaneously,  using  some  procedure  such  as  the 
two-dimensional  Newton-Raphson.  Further,  in  order  to  obtain  an  MLB  in  the  special 
case  when  there  are  no  censored  data,  merely  set  k «n.  Then,  from  equations  (3.4-1)  and 
(3.4-2)  and  by  setting  d(x)  « P(x)  • R(x), 


and 


need  simultanoous  solution  ot  obtain  MLE’s  for  f*.  and  y. 

It  can  be  shown,  for  the  case  of  complete  uncensored  samples,  that  the  statistics 

A A 

X and  V ■ ^ 

7 
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I? 


have  a joint  distribution  which  is  independent  of  the  parameters  n,y. 


By  definition  fi,y  satisfy  the  equations 


A A 


71  <2  n 

<f) 


Now  by  adding  and  subtracting  fi  and  letting  yi « (x{  • fi)/y  which  is  stochastically 
independent  of  in  distribution  yields  the  two  equations  in  T,V  namely; 

■ ■ + t)  . 0, 1 4 E (j;*  + T)d(!;‘  + t) 


Clearly  the  distribution  of  T,V  does  not  depend  upon  n,y  and  neither  does 


■‘t'J 


This  same  result  holds  in  the  case  when  a sample  is  truncated  at  the  kth  ordered 
observation. 


3.4.1  TES'rS  AND  ESTIMATION  B'ROM  GROUPS  OF  DATA 

Throughout  the  following  discussion  suppose  that  there  are  several  samples  hrom 
different  populations;  i.e.,  there  are  m such  samples, 


Xij,...,Xnj  lbrj“],..„m 


where  nj  observations  are  taken  fVom  the  jth  population  with  parameters  (ojiyji/ij)  for 
j » 1, . . m.  Now  compute  the  two  sample  statistics 


Xj  ® Z*  Xjj 

J »j  i-1  « 


(3.4.1-1) 


■i’  ‘Hpr  E 

J i“l 


(3.4.1 -2) 
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where  Xy  as  a normal  distribution  transformed  by  with  parameters 

(o^,7j,Pj)  fori=  1, . . nj 


then 


1 


(3.4.1-3) 


are  independent  identically  distributed  standard  normal  variates.  Since  is  odd  it 
follows 

and.ao 

• BXy=Mj.. _ 

and  thus  Xj  is  an  unbiased  estimate  of  /nj.  But  note  also  tl\ai 
By  defining 

Oj*»var(Xjj)  (3 ,4. 1-4) 


then 


aj*»7j"EW"(05V 


(3 .4. 1-5) 


and  it  is  known  that  an  unbiased  estimate  of  o-j®*  is  the  sample  standard  deviation 
squared  sj^. 

Thus  it  is  seen  that  this  family  of  distributions  separates  the  variance  into  the  product 
of  two  factors.  The  first  factor  is  the  scale  parameter,  y,  while  the  second  factor,  say 
B(a),  is  determined  by  the  flexure  parameter,  a,  and  the  family  characterization 
function  (u;  i.e.. 


OO 

B(a)  = Eto"' (ca.y  = J w“‘  (y)’d<l>(^) 

«OQ 
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Thus  the  equation^  valid  for  only  large  samples,  of 


forces  the  utilization  of  an  independent  method  for  estimating  aj. 
In  one  case,  namely  when 

uj(x) = sinh"‘  (x)  = 8n(x  + J x*+D 


B(j(«)  may  be  computed  explicitly. 
Let  X ~ N (0,«*),  recalling  that 


E(c^)  ® e 


then 


B^(a)  “ E(slnh  X) 


E(e»^-2+e-*^ 




and  finally 

BJa)  = ^(e««*-l) 


from  above 


thus  for  samples  large  enough  that  - o-j* 


1 

I 
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8.4.2  DISCUSSION  OF  EXPLORATORY  APPLICATION  AND  RESULTS 


Four  data  aamples  were  selected  for  initial  evaluation  of  the  ML  estimator  just 
described.  These  trial  cases,  all  of  size  n » 10  observations,  were  chosen  because  of  the 
dissimilarities  in  their  estimated  parameters  based  on  the  STAR  estimator.  These  data 
and  their  estimated  parameters  are  given  in  table  2;  the  behavior  of  the  estimator  for 
tire  first  data  sample  is  shown  in  figure  10.  It  is  interesting  to  compare  figures  7 and  10. 
In  both  cases  the  generated  values  of  the  partial  derivatives  dL/d/u  « 0 and  dUdy  « 0 
have  been  plotted  for  the  identical  fatigue  data  sample.  It  is  readily  apparent  that  in 
figure  7,  when  the  flexure  parameter  was  treated  as  an  unknown,  no  clear  cut  solution 
appears  to  exist  as  both  the  plotted  curves  run  more  or  less  parallel  with  each  other 
toward  some  asympototic  value.  However,  in  the  case  of  the  known  flexure  parameter 
a ■■  0.64,  there  is  an  obvious  singularity  resulting  in  unique  solutions  for  both  the 
location  (/x)  and  scale  (y)  parameters. 

Next  a censored  sample  comprising  16  failure  observations  and  32  censored  times  was 
investigated.  The  results  of  this  test  and  the  data  are  given  in  figure  11.  Again  a clear 
singularity  is  shown,  demonstrating  that  the  MLE  procedure  is  functioning  for  the 
censored  case  also.  It  was  also  noted  that  convergonce  was  rapid. 

The  ready  convergence  of  the  MLE  solutions,  when  given  the  value  of  the  flexure 
parameter,  has  been  convincingly  demonstrated.  Therefore,  the  data  previously 
investigated  to  establish  this  parameter  were  reevaluated  using  the  value  of  a » 0.64 
and  MLE*s  obtained  for  the  scale  and  location  parametors.  These  results  have  been 
tabulated  in  tables  3 and  4. 

A study  of  tables  3 and  4 shows  that  the  flexure  parameters  remain  almost  constant 
regardless  of  comparisons  made  between  life  lengthfi,  types  of  specimens,  or  loading 
procedures.  Considering  first  the  life  length  comparison,  it  is  seen  from  these  tables 
that  although  the  data  have  been  subdivided  into  groups  depending  on  the  magnitude  of 
their  respective  location  parameters,  the  resultant  weighted  average  flexure  value  for 
each  subset  is  almost  identical.  In  table  3,  the  weighted  average  8n-normal  standard 
deviations  for  tho  same  data  groups  are  also  shown  for  comparison,  and  it  is  seen  that 
the  trend  of  the  deviation  values  is  to  increase  with  increasing  values  of  the  location 
parameter.  In  table  4,  the  flexure  parameter  estimates  of  the  investigated  data  are 
grouped  according  to  specimen  types  and  loading  procedures;  once  more  it  is  noted  that 
this  parameter  remains  uninfluenced  by  such  variables.  However,  the  scale  parameter 
estimates  are  seen  to  have  more  variation  both  overall  and  when  comparing  specimen 
types  or  loading  procedures. 

These  observations,  therefore,  have  given  some  substance  to  the  initial  hypotheses  of: 
U)  a flexure  parameter  that  is  constant  and  independent  of  life  length,  geometry,  and 
loading  considerations;  and  (2)  a scale  parameter  that  reflects  geometry,  loading,  or 
other  such  variables.  It  must  be  pointed  out  that  this  study  has  not  included,  to  date, 
other  materials  such  as  titanium  or  steel  because  of  the  difficulty  in  readily  obtaining 
sufficient  suitable  data.  The  comments  just  noted  have  been  based  solely  on  aluminum 
alloy  information  and  may,  therefore,  require  some  adjustment  when  fiirther  data  are 
evaluated.  It  should  be  emphasized  here  that  aluminum  fatigue  data  with  both  very 
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little  and  extremely  large  scatter  were  included  in  this  investigation,  and  no  systematic 
difTerences  in  their  estimated  flexure  parameters  were  observed  when  compared  with 
estimates  from  more  typical  looking  test  results.  Therefore,  it  is  believed  that  the  value 
of  the  flexure  parameter  will  not  alter  by  much  for  materials  such  as  titanium  and 
steel. 

Tables  3 and  4 also  give  the  comparison  of  results  obtained  from  the  STAR  and  MLE 
estimat.jrs;  it  is  quite  clear  that  both  methods  produced  almost  identical  estimates  with 
regard  to  the  location  parameter,  n,  values.  In  the  case  of  the  scale  parameter,  y,  the 
difference  in  the  MLE  and  STAR  results,  although  still  small,  was  more  noticeable.  The 
scale  parameter  is,  to  some  extent,  dependent  on  the  magnitude  of  the  flexure 
parameter.  As  a result,  when  the  constant  MLE  flexure  value  is  less  than  or  larger  than 
the  STAR  estimated  flexure  parameter,  then  the  MLE  scale  parameter  is  larger  than  or 
lesser  than,  respectively,  the  STAR-estimated  scale  values  (see  table  2). 

Finally,  it  was  noted  that  there  was  a tendency  toward  a small  lowering  of  the  MLE 
scale  parameters  with  regard  to  the  STAR  scale  estimates.  This  can  be  seen  in  the 
results  given  in  both  tables  3 and  4.  Nevertheless,  the  overall  trend  of  the  results 
obtained  by  the  MLE  method  is  the  same  as  that  given  by  the  STAR  estimator,  and  the 
imposition  of  a fixed  value  for  the  flexure  parameter  does  not  appear  to  alter  the 
behavior  of  the  other  two  parameters. 

3.5  SIMULATION  METHODOLOGY  TO  GENERATE  SPECIFIED  r ORDERED 
OBSERVATIONS  OUT  OF  A SAMPLE  OF  SIZE  n 

The  structural  performance  of  a fleet  of  airplanes  is  typically  based  on  the  performance 
of  the  weakest  airframes,  with  the  result  that  the  extreme,  isolated  incidents  becoht^ 
the  basis  of  the  Information  for  subsequent  operation  and/or  modification  of  tho 
airplane.  Therefore,  consideration  must  be  given  to  the  important  problem  of  simulating 
fleet  performance;  i.e.,  the  generation  of  the  first  (early)  failures  (damage)  among  a 
larger  number  of  identical  components  that  are  expoiiod  to  the  same  operational  and 
loading  environments  and  that  are  assumed  to  be  independent.' 

One  way  to  do  this,  if,  for  example,  the  first  r observations  in  a sample  of  size  n were 
desired,  would  be  to  generate  the  n observations,  sort  them,  and  then  discard  all  but 
the  first  r ordered  observations.  It  can  be  seen  that  this  procedure  is  both  inelegant 
and  inefficient,  and  when  repeated  a large  number  of  times— typical  of  simulation 
procedures — it  becomes  very  costly.  Alternatively,  the  problem  can  be  formulated 
mathematically  as  follows. 

Generate  the  first  r ordered  observations  out  of  n,  given  the  density  f(x)  for  0 < x < m 
of  each  independent  observation.  Tlio  joint  density  of  the  first  r ordered  obsorvatiuns, 
say  (X|,  X2. . . .,  Xr),  is 

( n-r 
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where  the  distribution  and  survival  distribution  are  defined,  respectively,  by 

X 

fXx)=  f f(t)dt,  F(x)  = l-P(x),  fofx>0 
o 

By  setting  r - 1,  it  can  be  seen  that  the  marginal  density  of  the  first,  smallest 
observation  is 


i. 


g(x,)  = nf<x,)(P(x,)l 


n-1 


0<x,  < <» 


(3.6-2) 


'i 

I'r 


It  is  found  by  use  of  the  calculus  of  probabilities  that  the  conditional  density  of  the 
(j+l)st  ordered  observation,  given  the  value  of  the  jth,  is 

(nj)fl[Xj+i)[F(Xj4. ^ ^ (3.6-3) 


g(Xj+l)(Xj) 


checking  that 


[F()^)] 


n-i 


< Xj+1  < “ 


^I 


«(x,,x,,...,Xr)-g(x,)ri._j  J!(Xj+i|j^) 


which  is  in  agreement  with  equation  (3.6-1). 

Now  define  the  hazard  function  and  the  hazard  rate,  respectively,  by 


Il  = -CnP,h  = H' 


for  the  variate  with  survival  distribution  F. 

The  formula  for  the  generation  of  the  first  r ordered  observations  out  of  n,  with  r « n, 
in  terms  of  transformations  by  the  hazard  function  of  linear  combinations  of 
independent,  identically  distributed  exponential  variates  with  unit  mean,  is  now  given; 


If  H is  strictly  monotone,  and  so  has  an  inverse,  and  Y|,  Y^,  . . .,  Y|.  are  independent, 
identically  distributed  exponential  variates  with  unit  mean,  then 


Xk  = H 


(3.6-4) 


are  the  first  r ordered  observations  out  of  n independent  life  lengths,  each  with 
survival  distribution 


F(x)  = exp  I -H(x)|  for  X > 0 

I 
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By  letting  k = 1,  and  showing  that  Yi  ^ n H(Xi),  which  implies  that  Xj  has  the  density 
given  in  equation  (3.6-2),  it  can  be  found  that 

PpC,  >X,]a;<jxp  |-nH(x,)| 


Differentiating  this  to  obtain  the  density  of  X],  it  is  seen  that  the  density  is  that  given 
in  equation  (3.6-2). 


Supposing  that  equation  (3.6-4)  has  been  verified  for  k 1,2 j and 

4i 


(3.6-6) 


it  follov/B  fVom  the  conditional  density  of  equation  (3.6-3)  that  for  a given  xj 

Yj+i  = (n-j)[H(Xj.n)-H(x,)] 


in  exponential  with  unit  mean.  Thus, 


H (Xj+i)  = H(3^) 


and,  therefore,  from  equation  (3.6-6) 


Vi 

■ mil  ..  ..i.a 

n+l-i  n-J 


where 


IL. 

n+l-i 


is  regarded  us  being  given.  Therefore,  unconditionally 


It  is  worth  repeating  that  the  practical  method  for  eimulating  censored  data,  developed 
above,  applies  only  to  the  technique  of  censoring  by  number;  i.e.,  where,  after  the  first 
few  failure  observations  have  been  obtained,  the  data  gathering  is  suspended  or 
aborted.  Thus,  all  the  nonfailure  values  are  identical  and,  furthermore,  must  be  equal 
to  or  exceed  the  last  observed  failure  value. 

8.5.1  APPLICATION 

The  simulation  method  described  in  the  preceding  paragraphs  can  also  be  applied  to 
distributions  other  than  the  ones  being  used  in  the  following  examples.  For  consistency, 
it  has  been  assumed  here  that  X is  a {^-normal  variate.  Then,  by  definition 

P[X  > xl  = [- (^)]  - e-iiM  for  X > 0 


where 


• OO 


- eo  < X < «« 


is  the  standard  normal  distribution. 

Thus,  it  can  be  found  that 

H(x)  = -l>n«l>[-^{(^)]  x>0 
and 

H"’  (y)  = (e"’' )]  y > 0 


The  particular  case  of  interest,  as  before,  is  the  f-norma!  distribution  defined  as 

i(x)  = vf?r-(l/>/lc)  forx>0 

Consider  now  the  three  statistics,  T,  U,  and  V,  defined  as  follows: 
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which  are  known  to  be  invariant  with  regard  to  their  pararaotera.  The  diatributions  of 
these  atatistica  are  obtained  empirically  by  first  generating  censored  data,  using  the 
developed  MLE  methodology  to  estimate  the  parameters  and  then  computing  T,  U,  and 
V.  Finally,  this  procedure  is  repeated  many  times,  and.  the  results  are  ordered  to  form 
the  empirical  distributions. 


3.S.2  DISCUSSION  OF  EXPLORATORY  APPLICATION  AND  RESULTS 

ieveral  samples  of  censored  data  (ref.  16)  are  listed  in  table  5.  These  particular  data 
have  been  selected  as  examples  because  they  are  fairly  large  samples  and  also  are  of  the 
properly  censored  type.  It  is  noted  from  table  5 that  the  extent  of  censoring  ranges  from 
20%  to  90%,  and  that  sample  sizes  range  fVom  46  to  64  specimens. 

Figures  12  through  14  illustrate  some  results  obtained  by  using  the  simulation 
procedure  described  previously.  Figure  12  is  a simulation,  of  the  first  sample  listed  in 
table  5 and  presents  the  simulated  sampling  distributions  of  the  T and  V statistics  for 
a sample  of  size  52  where  79%  comprised  failure  observations  and  21%  were  nonfailure 
censored  values.  For  comparison,  the  empiric  distributions  of  the  T and  V statistics 
for  uncunsored  samples  of  the  same  size  as  the  number  of  failure  observations — in  this 
example  42  specimens — have  been  superimposed  in  the  figure  as  dotted  lines.  A study  of 
this  figure  shows  the  slight  increase  in  the  spread  of  the  solid  line  which  represents  the 
censored  case. 

Figures  13  and  14  present  simulated  distributions  representative  of  other  samples  from 
table  6.  A comparison  of  these  figures  clearly  illustrates  the  increase  in  the  amount  of 
bias  in  the  distributions  of  both  the  T and  V statistics  as  the  level  of  censoring 
becomes  larger.  This  bias  is  behoved  to  be  a result  of  Hist  the  extent  of  the  censoring  of 
the  data  and  second  the  limited  size  of  the  sample.  Again  returning  to  the  figures,  it  is 
readily  apparent  that  the  bias  in  each  uncensored  sample’s  V'Statistic  distribution  is 
less  than  that  exhibited  by  the  complementary  censored  sample.  Furthermort'v  it  can  be 
noted  that  bias  does  increase  as  the  sample  size  decreases,  except  for  the  T-statistic 
distributions  of  the  uncensored  cases  which,  as  expected,  appear  fairly  unbiased. 

Finally,  figure  15  is  presented  to  illustrate  the  relative  impact  on  the  T and  V 
distributions  of  the  absolute  number  of  failure  observations  and  the  level  of  censoring. 
This  has  been  achieved  by  maintaining  the  ratio  of  failure  to  suspended  observations 
but  doubling  the  sample  size.  For  this  example,  the  data  already  shown  in  figure  14 
were  selected  as  the  reference;  i.e.,  a sample  comprising  eight  failure  observations  and 
40  suspensions.  Now,  by  doubling  this  to  16  faihi'*es  and  80  suspensions,  the  level  of 
censoring  has  been  maintained  at  th&  original  83%,  but  the  absolute  number  of  failure 
observations  has  been  increusud  by  eight  data  points.  A study  of  figure  16  shows  that 
the  level  of  cemu  ttig  becomus  correspondingly  less  important  as  the  number  of  observed 
valueu  increases. 

One  additional  check  was  conducted  and  io  presented  in  figure  16.  It  should  be  noted 
that  all  the  preceding  ilguves  have  been  based  on  a flexure  parameter  (a  - 0.64). 
Therefore,  in  order  vo  check  whether  variation  in  the  value  of  this  parameter  should  be 
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a consideration,  the  first  two  data  samples  of  table  5 were  again  simulated  using 
different  values  of  a;  the  results  were  plotted  in  Hguie  16.  It  is  readily  apparent  from 
this  figure  that  significant  differences  in  the  distributions  were  not  observed. 

A direct  result  of  this  developed  and  demonstrated  simulation  methodology  is  the 
ability  to  obtain  bounds  for  parameter  estimates  made  from  censored  data.  This  is 
readily  apparent  from  the  figures  where  the  deviation  of  the  estimated  values  from  the 
true  values  is  clearly  illustrated.  For  example,  from  figure  12  can  be  seen  that  the  scale 
parameter  90%  bounds  are  0.76  at  the  5th  percentile  and  1.25  at  the  95th  percentile. 

3.6  DEVELOPMENT  OF  UNBIASED  MINIMUM  VARIANCE  ESTIMATORS 
FROM  THE  ML  ESTIMATORS  OBTAINED  PROM  CENSORED  FATIGUE  DATA 

The  following  discussion  is  based  on  the  assumptions  that  the  variate  X denotes  the 
logarithm  of  fatigue  life,  (i.e.,  X ~ Cn  N,  (see  sec.  2.1),  that  the  appropriate  transform 
(sec.  2.1.2)  is  taken  to  be 


oj“'  (X)  =Kn(x+Vx^  + l)  -«><x<oo 


and  that  the  shape  (or  flexure)  parameter  is  known,  but  the  locution  and  scale 
parameters,  and  y,  reopectively,  are  unknown  (sec.  3.3.1). 

Assume  that  a conjured  sample  of  the  first  r ordered  observations  out  of  n is  given. 
The  computational  procedures  described  in  section  3.4  are  known  to  solve  the  implicit 
ML  equations  necessary  to  obtain  the  MLE’s,  fi  and  Therefore,  fur  any  particular 
sample  of  size  (r,n),  the  statistics 


y 


V- 


T 

V 


are  known  to  be  invaiiant  with  respect  tc  the  parameters,  but  their  sampling 
distributions  must  be  found  by  the  simulation  technique  given  in  section  3.5.  Thus  by 
generating  a sequence  of  observations,  the  distribution  of  the  statistics 

(Tj,  Vj)  for ! = 1 , . . .,  ni 
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can  be  described  when  m is  large.  Furthermore,  the  sampling  bias  can  also  be  obtained 
by  computing  the  expectations,  valid  only  for  large  samples 


and 


m 


EV  = 


m 


i=I 


. 1 

& 


The  variance  of  these  statistics  is  given  by 

m 


H(T-abV)>  = i 2 ^ ^ TiV,  + ^ E VjV  - t* 


I I 

va,(V)4  E V|>.^=o> 

Thus,  for  any  real  sample  of  censored  fatigue  data  of  the  same  sample  size  (r,n)  as  the 
simulation,  the  ML  estimators,  ju  and  y can  be  used  to  construct  unbiased  estimators 
of  the  true  parameters  (j.  and  y by; 

]u  = 0 + ab^ 

^ “ b^ 


It  follows  that  these  ostimates  are  unbiased  since 

E (§  - M)  *=  7E  (^)  + at7EV  ==  0 
E$  «bTEV  =7 


3.6.1  DISCUSSION  OF  EXPLORATORY  APPLICATION  AND  RESULTS 

Seven  ’ censored  data  samples  from  reference  16  were  analyzed  using  the  techniques 
just  specified.  An  attempt  was  made  to  only  work  with  sample  sizes  that  were  both 
similar  and  large.  The  data  used  for  this  example  are  listed  in  table  5;  it  is  seen  that 
the  extent  of  censoring  varies  between  20%  and  90%.  Figure  17  presents  the  scale 
pernineter  sampling  bias  results  obtained  by  simulation  for  those  data  and  clearly 
illustrates  the  increased  bias  that  occurs  with  fewer  real  observations.  Two  extra  points 
have  been  spotted  on  this  curve  to  show  the  effects  of  decreased  total  sample  sizes.  In 
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each  case,  it  is  noted  that  the  bias  is  increased  with  the  smallest  sample  showing, 
relatively,  the  most  bias.  Figure  18  shows  the  location  parameter  sampling  bias  results 
for  the  same  set  of  data  samples,  and  the  trend  is  obviously  similar  to  that  noted  for  the 
scale  parameter. 

'fhese  initial  results  show  a promising  consistency  which  suggests  that  a matrix  of 
factors  might  well  be  a practical  goal.  This  would  then  circumvent  the  need  to  do  a data 
simulation  every  time  an  unbiased  estimate  is  made  of  the  distribution  parameters  from 
a data  sample. 


3.7  SIMULATION  OF  THE  DISTRIBUTION  OF 
ESTIMATED  POPULATION  PARAMETERS  OF  CENSORED  DATA 

Assume  that  a data  sample  of  size  (r,n)  is  one  which  has  been  truncated  at  the  r 
ordered  observation  out  of  n.  Given  values  of  r and  n,  a collection  of  m different 
data  sets  must  be  generated,  each  one  being  the  first  r out  of  n ordered  observations 
drawn  from  the  standard  distribution.  In  this  application  the  standard  distribution  is 
detined  by  g = 0 and  y - 1 . 

The  jth  sot  of  these  i values  is  denoted  by 

y,j Vq  forj  = l,...,  m 


Solving  for  the  MLE’s  of  the  location  and  scalo  parameters  for  these  values  allows  the 
determination  of  the  statistics 


for  j = 1 


Thus,  if  m is  iarge,  the  empiric  frequency  distributions  of  these  statistics  can  be 
described. 


Now  consider  a second  different  data  sample  of  size  (r^,n^ ).  By  repeating  the  above 
procedure  the  corresponding  statistics 

Tjt,Ujt,Vjt  forj-l m 


ar<5  computed. 

The  following  ratio 


Vj/Vjt  forj=l,...,  m 


is  calculated,  and  its  cumulative  frequency  distribution  is  determined. 
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Thus,  it  can  be  seen  that  by  this  process  the  joint  bounds  of  the  statistics  Vj  and  Vy 
have  been  generated.  Furthermore,  it  follows  that  if  the  two  censored  samples  of  size 
(r,n)  and  both  have  the  same  true  scale  parameter  (i.e,,  y then  their  ratio 

is  simply  the  ratio  of  estimated  scale  parameters  y and  (y)'!*.  Therefore,  by  calculating 
the  actual  ratio  of  the  scale  parameter  estimaUs  for  the  two  samples  and  comparing 
this  ratio  with  the  simulated  distribution,  it  is  possible  to  test  the  hypothesis  that  tho 
two  sets  of  data  come  from  a population  with  the  same  scale  parameter  (i.e.,  with  the 
same  unknown  y,  but  possibly  different  m’b). 

For  example,  assume  that  it  is  necessary  to  test  the  hypotheses  of  equivalent  scale 
parameters  at  the  10%  level;  i.e.,  to  test  whether 

Hq:  7 = 


Now  compare  the  ratio  of  the  MLE*s,  y and  y*^,  for  the  two  real  data  samples  with  the 
simulated  distribution,  and  note  whether  the  ratio  falls  outside  the  5th  and  96th 
percentiles.  If  it  does.  Ho  is  rejected  at  the  10%  level,  asserting  a 90%  certainty  that  the 
truo-acale  parameters  of  the  two  data  samples  are  different. 

A similar  process  to  test  the  hypothesis  that  the  two  samples  (r,n)  and  (r'^,n‘l‘)  also  have 
the  same  true  location  parameter  /a  has  been  developed  for  application  to  those 
samples  that  pass  the  y test.  Here,  a simple  relationship  between  the  MLE’s  of  the 
samples*  location  parameters  will  not  suffice  as  the  prior  test  has  already  demonstrated 
that  both  samples  most  likely  have  the  same  true,  but  unknown,  scale  parameter  y. 
Thus,  a comparison  should  be  made  between  estimates  of  the  two  samples’  location 
parameters,  which  have  been  based  on  this  true  scale  parameter.  Consequently,  the 
empiric  distribution  of  the  "T”  statistic  (T  =•  (n-flVy)  is  of  interest  for  the  location 
parameter  hypothesis  test.  The  following  relationship 

Tj-Tjt  forj  = 1, 2, . . .,  m 


is  calculated,  and  its  cumulative  frequency  distribution  is  determined  in  a similar 
manner  to  the  V ratio  already  described. 

The  value  of  the  scale  parameter,  which  must  be  common  to  both  data  sumples  (r,n)  and 
(r''',n‘l’)  in  order  to  test  the  hypothesis  of  equality  of  the  location  parameters,  is 
computed  in  the  following  manner.  The  empiric  distributions  resulting  from  the 
simulation  of  the  V and  statistics  are  used  to  obtain  the  bias  factors  b and  b*^ 
and  the  spread  factors  s and  sf*  which  are  the  spread  of  the  V and  distributions, 
> espectively,  at  some  fixed  percentile.  Recalling  the  presumption  that  y yt«y  which 
is  unknown,  the  minimum  variance  unbiased  estimate  of  y,  based  on  the  statistics  y and 
y+,  is  given  by: 

f + ct^t  = cb^  + etbt^t 
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where 


Cj  = (bjSj)“V  2 ^hjSj)  ^ i = I,  2 

j=:| 


where  Cj  c,  C2  = ct 


Note  that  e ia  aBsumed  to  be  proportional  to  «r  the  standard  deviation  of  the  V 
statistic. 

This  estimate  can,  of  course,  bo  generalized  for  the  case  yj  =>  yj  “ fo** 

K>  2.  Under  those  conditions 

1!  C 

Si = 7V 
I 1 


This  is,  in  fact,  the  minimum  variance  estimator  since  by  the  Cauchy  inequality;  i.o.. 


it  follows  that 


£ (OjhiSi)^  • 2]^  (hjSj)  -*  > I 


for  any  set  of  ci  0 so  that 


2 

Tf-r' 


It  should  bo  noted  that  the  spread  s urs*^  rather  than  the  variance  has  been  adopted 
for  the  sake  of  computational  expediency.  However,  it  is  presumed  that  the  minimum 
variance  unbiased  estimate,  y,  will  produce  a common  value  so  close  to  the  true  valdo  of 
y that  the  empiric  distribution  of  T cun  be  used  justinubly  in  this  investigation  for 
tests  or  confidence  bounds  on  the  location  parameter /a.  Although  this  approximation 
may  not  be  exact,  the  actual  difference  is  unlikely  to  vitiate  any  of  the  conclusions 
reached  while  improving  the  basic  capability  of  the  reliability  analysis  method  to  treat 
censored  data. 
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3.7.1  DISCUSSION  OF  EXPLORATORY  APPUCATION  AND  RESULTS 


Several  samples  of  censored  fatigue  data  from  table  6 have  been  paired  for  comparison 
according  to  the  methods  aescribed  in  the  previous  paragraphs,  and  the  results  are 
presented  in  figures  19  through  22. 

Figure  19a  and  19b  present  the  simulation-generated  cumulative  frequency 
distributions  of  the  V and  T joint  statistics  for  two  censored  samples  of  size  (41,62) 
and  (33,53).  From  these  figures  it  is  clear  how  bounds  for  the  "V-ratios"  and 
"T-differences”  can  be  established.  For  example,  it  can  be  seen  that  the  90%  bounds  for 
the  V-ratio  for  two  samples  of  size  (41,62)  and  (33,53)  are 


In  a similar  way  the  values  -0.26  and  0.232  are  the  90%  bounds  for  the  joint  statistic 

(T41,52-T33,53). 

Table  6 gives  the  MLE’s  of  the  scale  and  location  parameters  for  the  various  data 
samples  examined.  It  is  noted  that  the  estimated  scale  parameter  of  the  first  test 
sample  (of  size  41,52)  has  a value  of  » 0.40  and  the  second  test  sample  (of  size  33,53) 
has  a value  of  ^2  = 0.696.  Therefore,  the  ratio  (^i/^'>)  “ 0.676. 

Returning  to  the  results  of  figure  19,  discussed  earlier,  it  was  noted  that  the  low  value 
of  the  90%  bound  of  the  V-ratio  was  0.696.  Therefore,  it  can  be  concluded  that  it  is  90% 
certain  that  the  true  scale  parameters  of  the  censored  data  samples  1 and  2 are 
different. 

Referring  once  more  to  table  6,  it  is  seen  that  sample  3 is  also  the  same  size  as  sample  2 
(i.o,,  r“33,  ns'eS),  and  so  the  simulation  results  are  equally  applicable  to  either  the 
(sample  1/sample  2)  or  (sample  1/sample  3)  ratios.  From  table  6 it  is  found  that 
^sample  a = 0.636  and,  therefore,  the  ratio  (^1^3)  “ 0.748.  Thus,  in  this  case  the  teat 
data  ratios  fall  within  the  90%  tolerance  bounds  of  figure  19,  and  the  hypothesis  that 
the  uame  truo-scale  parameter  applies  to  both  samples  1 and  3 cannot  be  rejected  at  this 
level. 

Figures  20  and  21  illustrate  the  joint-scale  and  joint-location  statistics  obtained  through 
simulation  of  two  other  combinations  of  samples.  Figure  20  is  interesting  because  it 
represents  the  expected  variability  between  two  independent  and  identically  sized 
samples,  such  as  items  2 and  3 of  table  6,  From  figure  20,  it  cun  be  found  that  the 
hypcthesis  of  the  scale  equality  of  samples  2 and  3 cannot  be  rejected  at  the  90%  level. 

Figure  21  applies  to  data  samples  4 and  6,  and  from  those  results  the  hypothesis  of 
scale  equality  of  samples  4 and  6 is  not  rejected  at  the  previously  selected  test  level  of 
90%.  This  figure  is  interesting  because  sample  6 ia  over  90%  censored.  Consequently, 
there  is  considerable  bias  to  the  estimates  of  Its  scale  and  location  parameters  as  is 
apparent  ir  Uio  "V”  and  "T”  plots  of  figure  21. 
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This  discussion  has  covered  the  comparison  of  paired  samples  based  on  the  estimated 
scale  parameters,  and  for  the  examples  tested,  it  was  found  that  the  sample  pairings 
such  as.  1 and  3,  2 and  3,  and  4 and  3 were  not  rejected.  The  minimum  variance 
unbiased  estimate  of  the  scale  parameter  of  each  of  these  three  paired  samples  was 
calculated,  and  then  the  MLE  of  the  location  parameters  of  each  sample,  with  specified 
flexure  and  scale  parameters,  was  computed.  These  results  are  shown  in  table? 
togother  with  the  diiTerences  of  the  revised  location  parameters.  The  difference  values 
were  then  compared  to  the  location  parameter  bounds  shown  in  figures  19,  20,  and  21;  it 
was  clear  that  the  hypothesis  of  location  equality  for  each  paired  sample  could  not  be 
rejected  at  the  90'Jt>  level  selected  earlier. 

It  is  believed  that  the  few  examples  discussed  in  the  preceding  paragraphs  are  sufllcient 
for  illustrating  the  developed  procedure  for  testing  whether  or  not  censored  samples 
could  have  come  from  the  same  parent  population.  It  has  been  noted  that  the  testing 
procedure  is  able  to  either  reject  or  not  reject  the  hypothesis  of  scale  equality  and 
subsequent  location  equality. 


4.0  SIGNIFICANCE  OF  MIXED  POPULATIONS 


4.1  THE  INITIAL  FLAW  CONCEPT 

The  structural  deaign  goals  and  subsequent  operation  of  the  airplan^i  are  based  on  the 
rationale  that  each  structure  is  manufactured  similarly.  The  performance  of  each 
structure  should,  therefore,  be  consistent.  However,  the  tolerances  that  are  necesaary 
both  during  the  production  of  the  material  and  its  subsequent  nianiiiiacture  into  n 
structural  detail  or  component  do  cause  a variation  in  the  performance  of  individual 
structures.  This  variability  is  usually  taken  into  account  during  the  design  by  the 
assumption  of  some  average  level  of  performance  which  is  then  downgraded  through  use 
of  various  factors  to  a more  represontative  performance  level. 

Therefore,  in  a roundabout  manner,  procedures  such  as  the  above  presuppose  the 
existence  of  a statistical  population  which  can  be  simply  described  by  central  tendency 
and  variability  parameters.  Tlius,  for  example,  the  fatigue  performance  of  any  structure 
reflects  the  performance  in  terms  of  time  or  cycles  to  o dotectabio  crack  of  the  typical 
structure,  i.e.,  it  comes  fVom  the  same  parent  population. 

However,  there  arc  data  available,  such  as  that  collocted,  analyzed,  and  presented  in 
reference  17,  leading  to  the  hypothesis  that  a numhor  of  early  structural  mcldenta 
belong  to  a population  that  is  both  different  and  worse  than  the  population  from  which 
the  typical  structures  come.  This  means  that  if  sufficient  data  were  available,  the 
frequency  of  plotted  times-to^crack  initiation  would  appear  bimodal  with  the  lower  mode 
representing  the  fatigue  observations  resulting  from  initial  flaws,  manufacturing 
damage,  etc.,  and  the  second  mode  representing  the  average  fatigue  performance  of  the 
typical  structure.  It  follows  then  that  the  present  inetho<l  of  fatigue  performance 
prediction,  based  on  the  behavior  of  the  typical  structure,  cannot  predict  the  very  early 
incidents  that  occur  as  the  result  of  an  inadvertent  flaw.  Under  this  circumstanca.  the 
analysis  of  the  reliability  of  a structure,  based  on  the  single  typical  population,  will  be 
unconaervativo,  especially  during  the  early  life  of  the  structure  prior  to  a routine 
inspection.  Therefore,  it  Is  important  to  recognize  the  initial  flaw  population  so  that  the 
inspection  procedures  can  be  adijusted  accordingly  and  a high  level  of  reliability 
maintained.  Furthermore,  the  current  Air  Force  philosophy  (ref.  16)  is  oriented  to  the 
concept  of  the  initial  or  preflaw  for  the  analysis  of  damage-tolerant  structure.  However, 
not  all  details  are  necessarily  proflawod,  and  one  of  the  tasks  of  this  study  is  the 
extension  of  the  reliability  analysis  methodology  to  assess  the  effect  of  the  presence  of 
details  with  a reduced  level  of  fatigue  performance  in  a fleet  of  aircraft.  Hence,  for  the 
purposes  of  this  investigation,  the  possibility  is  considered  for  some  structural  details  or 
components  entering  service  with  an  existing  flaw  or  damage  of  magnitude  equal  to  or 
greater  than  some  specified  minimum  detectable  size. 

The  model  that  has  been  developed  accounts  for  both  the  probable  frequency  of  flawed 
structures  entering  service  and  the  probable  flaw  or  damage  size  in  the  structure.  It  is 
realized  that  there  may  be  considerable  difficulty  in  obtaining  suitable  information  on 
the  likelihood  of  flaws  or  damage  in  new  structure  and  the  probable  magnitude  of  these 
flaws  or  damage. 


Nevertheless,  work  such  as  that  of  reference  17  already  provides  indications  that 
applicable  input  information  can  be  obtained.  If  not,  perhaps  judicious  selection  of  input 
values  m^y  be  done.  Consequently,  the  development  of  an  intial  flaw  model  and  its 
inclusion  in  a structural  reliability  analysis  system  should  provide  an  added  measure  of 
realism  to  the  ensuing  results. 
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5.0  THE  RELIABILITY  MODEL 


5.1  STRUCTURAL  RELIABILITY  ANALYSIS  METHOD 

An  analysis  method  to  estimate  the  structural  reliability  of  an  airplane  structure  based 
on  the  interaction  of  cumulative  and  maximum  operational  loads,  residual  strength  of 
fatigue  damaged  structure,  and  routine  inspection  and  crack,  detection  and  repair  of  the 
structure  was  developed  in  earlier  investigations  and  reported  in  references  12  and  13. 
Some  of  the  design  variables  that  were  included  in  this  reliability  system  were  the 
central  tendency  time,  to  initiation  of  a detectable  crack  and  its  var.' ability,  the  size  of 
the  detectable  crack,  hnd  the  central  tendency  crack  growth  rate  and  its  variability. 
Monte  Carlo  simulation  techniques  are  then  used  to  generr.te  a timo  to  a detectable 
crack  of  specified  size  and,  subsequently,  to  grow  the  crack  until  it  is  either  arrested  by 
a crack  stopper  or  the  part  is  failed. 

This  reliability  system  was  based  on  the  principle  of  the  .lingle  typical  population,  as 
discussed  in  section  4.0,  from  which  random  initiation  times  are  selected.  Obviously, 
such  a model  could  guiieratu  some  crack  initiation  times  approaching  or  at  "zero  time" 
depending  on  the  values  of  the  central  tendency  and  variability  parameters  initially 
input  into  the  reliability  system.  However,  this  would  be  as  a direct  result  of  the  typical 
component  population’s  tail  and  not  representetlvo  of  the  nontypical  structure  with  the 
flaw.  Consequently,  it  was  decided  to  expand  the  reliability  system  of  references  12 
and  13,  and  thereby  increase  'its  comprehensiveness  to  include  the  initial  or  preflaw 
model  discussod  in  the  following  pages. 

FLAW  SIZE  DENSITY 

Given  that  flaws  of  differing  sizes  are  distributed  randomly  in  :.,tructural  materials,  it  is 
assumed  that  this  flaw  size  density  is  exponentially  decreasing. 

Lot  T denote  Iho  random  initial  flaw  size  with  parameters  p.^  and  t,  related  by  the 
constraint 

l’[f  < ’•»]  -•  fr 


Here  to  is  some  practical  minimum  detectable  flaw/damugu  size,  and  p^,  0 <p^  < 1,  is 
that  portion  of  the  area  under  the  density  curve  representutlve  of  the  flaws  that  are 
below  T()  in  size  and  therefore  not  detectable. 


Now  assume: 
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where  Z ie  the  standardized  normal  variate  with  zero  mean  and  unit  deviation  and  a is 
the  correlation  parameter  between  the  minimum  detectable  flow/damage  size  and  its 
frequency,  as  represented  by  the  following  relationship: 

a = T„/f/ 

where  (o^is  the  standard  normal  quantile  at  a given  probability  level,  a,  with  a taken  as 
(l+p^/2;  and,  v is  the  magnitude  at  which  flaw  size  density  in  the  structure  is 
exponentially  decreaaing. 

Flaws  or  damage,  in  a structural  component,  of  size  <ru  are  considered  to  be  typical, 
and  the  standard  operational  procedure  is  to  compute  the  time  to  initiation  of  a crack  of 
size  TO  and,  subsequently,  to  determine  crack  growth  to  critical  size.  However,  when 
flaws  or  damage  of  size  *rQ  exist,  tho  structural  component  is  considered  to  be 
preflawed  and  only  the  subsequent  crack  growth  in  the  structure  need  be  considered. 
The  procedure  for  tills  oporatio'i  follows: 


Generate  Z - N(0,1);  test  whether  Z < where  a ia  the  probability  level  as  before. 


Now  if  this  is  the  case,  then  no  flaw  or  damage  of  detectable  size  exists  in  the  structure 
and  the  time-to-crack  initiation  must  be  computed.  However,  if  Z > then  a flaw  or 
damage  of  size 


is  already  present  in  tho  structural  component  so  that  only  the  subsequent  crack  growth 
in  the  structure  need«  to  be  oonsldorod. 


A 


From  this  procedure,  it  followc 


PIT  < 


and  a continuous  density  for  flaw  sizes  exceeding  ro  of: 


A 1 o-l/J 

=V  J 1 2. for  r > r„  \ i>  > 0 


Now,  let  n be  the  size  of  a largo  initial  flaw  or  crack  which  has  some  small  probability 
e of  existing  in  a component  entering  service.  This  is  given  by  the  following  constraint 
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or 

PIa|Z|>'  = r,l  = e 


thus 
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Thin  reduces  to 
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where; 


and  solving  for  v gives 


2n(r./7;,) 


The  behavior  of  the  density  of  the  initial  flaw  sizes  exceeding  tq  varies  radically  with 
the  value  of  v > 0,  as  shown  in  the  following  schematic; 
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Proof:  In  all  cases,  for  given  v > 0 

F(r:»^)  = F[T  < rj  -P[(Z|  < --2«^[(^•A-l 


where  4>(x)  is  the  cumulative  distvibutior  of  the  standard  normal.  It  is  checked  directly 
that  F(to:i>)  = p,j.  and  by  differentiation  that 

T_  1 - IIP'  1 

I- 

which  reduces  to  the  preceding  result. 
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6.0  INTERACTION  OF  MATERIAL/ 
STRUCTURE  VARIABILITY  WITH  RELIABILITY  LEVEL 
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6.1  DISCUSSION  OF  EXPLORATORY  APPLICATION  RESULTS 

The  renults  of  parametric  studies  using  tho  expanded  structural  reliability  analysis 
system  are  discussed  in  this  section.  The  hypothetical  cases  that  have  been  investigated 
are  based  on  the  following  assumptions: 

• The  characteristic  life  to  a detectable  crack  of  size  0.03  in.  is  120  000  hr. 

• The  minimum  or  threshold  size  for  crack  detection  is  0.02  in. 

• The  scheduled  inspection  period  is  7500  hr. 

4 The  fail-safe  residual  strength  of  a component  with  a cracked  detail  is  taken  to  be 

limit  load  capability. 

• The  reference  stress  intensity  factor  for  aluminum  is  taken  as  16  ksi  vln. 

9 The  crack  growth  rate  for  aluminum  is  taken  as  1.25  x 10'^  in./oyole. 

• The  critical  stress  intensity  factor  for  aluminum  is  taken  as  75  ksi  Vim 

• The  reference  stress  intensity  factor  for  steel  is  taken  as  20  ksi  VTn". 

• Tho  crack  growth  rate  for  steel  is  taken  as  6.6  x 10  in./cycle. 

• The  critical  stress  intensity  factor  for  steel  is  taken  as  100  ksi  VinT 

It  should  be  noted  that  any  deviations  from  these  values  are  indicated  in  the 
appropriate  figures  (23  through  30). 

Two  different  loading  cases  wore  also  used  to  describe  the  mean  rates  of  occurrence  of 
specific  structural  load  levels  per  hour.  These  loud  cases  are  illustrated  in  figure  22  and 
represent  exceedance  values  that  may  be  expected  at  two  different  locations  of  a 
four-engined  jet  tanker/transport  type  airplane  on  a routine  mission. 

Figure  23  illustrates  the  results  of  comparisons  made  between  several  cases  with 
differing  preflaw  assumptions.  The  example’s  structure  was  assumed  to  be  basically  an 
aluminum  skin  panel  with  crack  stoppers  placed  at  14~in.  intervals.  The  preflaw  or  the 
initiated  fatigue  crack  was  assumed  to  bo  located  at  a central  hole  between  the  crack 
stoppers.  The  component  was  assumed  to  be  fail-safe,  so  that  if  a crack  propagated  to 
the  full  width  of  a bay,  the  remaining  structure  would  still  be  capable  of  carrying  limit 
load.  The  approach  used  to  describe  the  preflaw  density  parameters  was  determined  by 
selecting  two  flaw  sizes  and  their  respective  probability  levels  and  then  computing  the 
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exponent  that  would  best  fit  these  points.  The  reference  curve  at  the  lower  right  corner 
of  figure  23  represents  the  probability  of  failure,  or  alternatively,  the  reliability  of  the 
fail-safe  component  (just  described)  when  exposed  to  the  case  A loads  of  figure  22.  The 
other  curves  plotted  in  figure  23  illustrate  two  things;  that  component  failure 
probability  increases  as  the  probability  of  an  extant  detectable  flaw  at  "zero-time” 
increases;  and  also  that  the  probable  size  of  the  undetected  flaw  is  a signiilcunt 
parameter,  since  the  component  reliability  decreases  with  increased  flaw  size.  Another 
expected  observation  is  that  the  decreased  component  reliability  which  results  from  an 
assumption  or  a preflaw  (or  damage)  condition  is  most  noticeable  at  the  lower  lifetimes 
where  several  orders  of  magnitude  differentiate  between  the  flawed  and  unflawed  cases. 


Figure  24  illustrates  the  influence  of  the  residual  strength  on  the  failure  probability  of 
the  component.  Basically,  th^  two  lower  curves  represent  fail-safe  structures,  whereas 
the  upper  curve  is  for  a structural  component  whose  cracked  strength  is  little  more  than 
the  strength  of  the  detail  represented  by  the  single  bay.  In  these  examples,  the 
component  was  assumed  to  be  free  of  flaws.  However,  figure  2b  presents  the  same  cases 
except  that  a high  likelihood  of  preflaws  was  assumed.  It  is  clear  that  the  combination 
of  high  preflaw  probability  and  lov/  fail-safe  strength  results  in  an  unreliable  structure. 


Figure  26  is  included  here  to  illustrate  the  importance  of  the  load  environment.  The 
results  shown  in  this  figure  come  from  the  assumption  that  the  structure  is  exposed  to 
the  case  B loads  of  figure  22.  Consequently,  by  comparing  the  results  in  figure  23 
with  those  of  figure  26,  the  effects  of  the  load  cases  A and  B,  respectively,  on  structural 
reliability  can  be  noted,  it  is  clear  from  this  comparison  that  the  load  exposure  term  is 
an  important  one;  it  can  alter  not  only  the  probability  of  failure  but  also  the  trend  of 
the  results.  Fur  example,  in  contrast  to  the  results  of  figure  23,  figure  26  does  not  show 
the  ve^y  large  differences  in  structural  reliability  between  the  unflawed  and  preflawed 
examples. 

Figure  27  illustrates  the  results  of  similar  studies  to  those  of  figure  26  except  that  in 
these  cases  the  type  of  structure  is  changed.  For  these  examples,  the  structural 
component  comprises  several  pinned  details,  so  that  loss  through  cracking  of  any  one 
detail  still  loaves  the  component  with  sufficient  residual  strength  to  sustain  limit  load. 
It  should  be  noted  that  the  various  parameters  used  to  obtain  the  results  shown  were 
selected  specifically  to  maintain  a similar  probability  of  failure  curve  for  the  reference 
unflawed  condition  as  shown  in  figure  26  for  the  skin-type  structure.  A comparison  of 
these  two  figures  clearly  illuatrates  those  results. 

Figure  28  represents  cases  identical  to  those  of  figure  27  except  that  the  material  has 
been  changed  from  aluminum  to  steel.  Once  again,  in  order  to  have  a reference  baseline 
for  comparison  with  other  cases,  the  parameters  wore  selected  to  produce  a probability 
of  failure  curve  for  the  unflawed  component,  similar  to  those  of  figures  26  and  27.  The 
results  of  the  equivalent  aluminum  and  steel  cases  show  that  the  latter  structure  is 
somewhat  more  sensitive  to  the  preflaw  assumptions. 
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Figures  2D  and  30  compare  the  ett'ects  of  the  residual  strength  assumptions  on  the 
aluminum  and  steel  structural  components,  respectively.  Although  it  is  clear  that  the 
trend  of  both  figures  is  similar,  it  is  equally  obvious  that  the  latter  steel  structure  is 
rof'ta  uensitive  to  the  residual  strength  term.  However,  in  both  cased  it  is  also  clear  that 
maintaining  a high  level  of  residual  strength  in  the  precraeked  state,  by  incorporating 
fail-safe  design  procedures,  is  very  important. 

It  id  obvious  Orom  these  uxaihples  and  discussions  that  the  reliability  anlaysis  method 
c.'^n  consider  and  evaluate  the  influence  on  ruliability  of  possible  flaws  in  a liew 
structure,  Furthermore,  the  presence  of  initial  flaws  or  even  other  conditions  leading  to 
mixed  (o<^  heterogeneous)  populations  of  fatigue  performance  in  terms  of  initiation  time 
to  detectable  crack  has  a definite  effect  on  the  fatigue  behavior  and  reliability  of  a fleet 
of  aircraft.  The  combination  of  the  probabilities  of  an  initial  flaw  and  the  usually 
expected  flaw-free  or  damage-free  structural  detail  will  tend  to  develop  an  S-shaped 
character  to  the  cumulative  distribution  ilmction  of  a group  of  details  as  compared  to 
the  usual  linear  function  when  plotted  on  normal  probability  paper  with  logarithmic 
abscissa.  If  there  is  a positive  probability  of  initial  flaws  or  cracks  at  time  zero,  the 
distribution  function  would  become  asymptotic  to  this  probability  level,  On  the  other 
hand,  should  there  be  a positive  waiting  time  before  flaw  or  crack  initiation  occurs,  a 
. reversed  flexure  would  cause  the  initial  toe  of  tho  distribution  to  rise  more  steeply  than 
the  expected  linear  representation.  Certainly,  the  inclusion  of  the  preflaw  or 
mixed-population  concept  into  the  reliability  analysis  scheme  is  important  to  the 
development  of  inspection  control  plans  and  maintenance  schedules  of  an  aircraft 
structural  system.  As  another  related  point,  the  effects  on  structural  reliability  of 
cracks,  which  originate  either  from  u preflaw  or  initiate  in  unflawed  structure  due  to 
fatigue,  are  sensitive  to  a number  of  interdependent  variables  including  loads, 
materials,  structural  conflgurations,  residual  strength,  and  crack  detectability.  Finally, 
it  should  be  noted  that  with  the  inclusion  of  fail-safe  components,  the  structural 
reliability  of  a system  can  bo  maintained  at  a high  level  even  with  the  assumption  of  a 
high  probability  of  a preflaw  in  a detail. 
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7*0  STRUCTURAL  DATA  INPUT  CONSIDERATIONS 
FOR  APPLICATION  OF  RELIABILITY  MODEL 
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The  durability  and  integrity  of  a structural  ayetem  exposed  to  a fatigue  environment 
are  responsive  to  the  application  of  reliability  analysis  technology.  However,  the  degree 
of  responsiveness  of  the  analysis  is  really  dependent  upon  the  modeling  details  and  the 
input  data.  Furthermore,  almost  every  rational  scheme  which  quantifies  the  elements  of 
the  fatigue  process  provides  some  measure  of  insurance,  though  not  necessarily  security, 
against  the  serious  consequence  of  unexpected  fatigue  damage  The  more  realistic  the 
approach,  the  more  insight  provided  for  the  task  of  determining  a numerical  measure; 
i.e.,  the  probability  of  structural  survival,  to  judge  progress  toward  operational 
structural  goals.  While  this  remark  suggests  that  the  application  of  reliability  analysis 
technology  alone  may  not  be  suiTicient  to  resolve  the  fatigue  problem,  the  technique  has 
a potential  to  provide  a workable  and  useful  tool  for  aiding  in  the  resolution  of  the  total 
fatigue  assessment  process  from  design  through  uaofVil  life. 
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Considerations  leading  to  this  somewhat  paradoxical  comment  are  those  related  to  the 
complexities  of  the  fatigue  process  itself,  as  well  as  the  applied  reliability  technology. 
Primarily,  cumulative  fatigue  damage,  whether  dealing  with  crack  initiation  or  stable 
crack  growth,  the  residual  strength  of  the  fatigue-damaged  structure,  and  the  coincident 
actual  mechanical  and  physical  environment  have  shared,  to  various  degrees,  in 
confounding  the  durability  and  integrity  goals  of  fatigue  critical  structures.  During  the 
development  of  the  reliability  analysis  scheme  of  this  study  and  the  preceding 
exploratory  work,  the  emphasis  has  been  toward  relating  the  reliability  analysis 
methods  to  recognized  or  measurable  structural  fatigue  properties  of  materials  in  their 
structural  form.  B'or  instance,  different  structural  details  of  the  same  material  can  have 
widely  differing  fatigue  performance  levels.  Furthermore,  idontical  details  in  differing 
loading  environments  also  can  exhibit  a different  response.  Recognition  and  separation 
of  these  complex  elements  are  believed  necessary  to  make  the  application  of  reliability 
technology  u tractable  tauk.  Skill,  understanding,  experience,  and  even  prescience  are 
some  of  the  ingredients  that  influence  the  degree  of  success  in  making  the  cumulative 
damage  prediction  process  match  operational  results.  Heneo,  this  reliability  analysis 
plan  is  directed  toward  utilizing  the  known  and  reasonably  identifiable 
matorial/structures  fatigue  properties.  Although  variability  of  Identical  structures 
under  a Hpoclfic  load  environment  is  universally  observed  to  respond  to  statistical 
treatment,  the  exact  nature  of  that  statistical  behavior  is  not  nearly  as  well  defined, 
particularly  in  reference  to  extreme  behavior.  Novertheloss,  this  approach  offers  a 
moans  to  explain  an  observable  behavior  (viz,  the  initiation  of  fatigue  damage)  in  a 
rational  manner  through  reliability  analysis  techniques.  As  the  important  technologies 
and  knowledge  pertinent  to  the  cumulative  damage  process  are  extended,  such  a 
reliability  analysis  scheme  will  incorporate  those  advances, 


As  an  example  of  the  preceding  concern,  the  static  strengths  of  commercially  available 
materials  are  treated  by  reliability  methods.  First,  variability  is  recognized  and  then 
characterized  by  presuming  a normal  distribution  for  the  specific  structural  properties. 
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such  as  ultimate  tensile  strength  or  tensile  yield  stress.  A confidence  level,  reflecting 
input. knowledge  (viz,  the  number  of  tests  performed  to  obtain  the  specific  value),  has 
been  established  for  certain  materials  such  as  the  aluminum  alloys.  With  that  baseline, 
a probabilistic  level  is  selected  to  determine  a strength  level  that  will  be  equaled  or 
exceeded  by  that  particular  alloy.  The  "A**  and  "B’'  material  property  values  of  the 
MIL-HDBK-6  compilation  of  such  structural  data  illustrate  this  application  of 
reliability  methodology.  The  universal  acceptance  and  demonstrated  usability  of  this 
technique  are  believed  reasonable  justification  for  the  procedure  and  encourage  its  use 
in  the  fatigue  evaluation  process.  However,  one  difficulty  in  this  phase  of  the 
application  of  reliability  methodology  to  fatigue  concerns  the  added  dimension  of  fatigue 
performance.  The  collection  of  the  static  properties  of  a material  is  a rather 
straight-forward  process.  For  instance,  10  static  test  specimens  from  each  of  10  different 
heats  of  a material  can  be  tested  at  a reasonable  cost.  The  mere  introduction  of  load 
cycling  and  the  very  local  nature  of  fatigue  damage,  initiation  in  structural  details, 
particularly  in  complex  components,  make  the  task  several  orders  of  magnitude  greater 
in  specimen  and  test  cost  as  compared  to  characterizing  the  static  strength  properties. 
However,  with  definition  of  the  nature  of  the  distribution  for  variability  of  the 
particular  fatigue  properties,  statistical  estimation  procedures  for  small  samples,  such 
as  one  or  two,  etc.,  provide  a means  to  estimate  the  tolerances  in  performance 
unavoidably  due  to  such  limited  Input  data. 

Continuing  the  static/fatigue  comparison,  the  determination  of  the  static  strength  of 
structures  from  reliability-determined  material  properties  introduces  a further 
complication.  Full-scale  static  tests  of  structures  made  from  known  materials  do  not 
always  fulfill  their  static  analysis  goals.  A natural  step  is  to  compare  the  results  of  this 
structural  testing  with  the  estimated  goal.  The  discrepancies  are  then  presumed  to  be 
variable  or  statistical  quantities  that  can  be  treated  in  a probabilistic  manner  and  then 
introduced  into  the  strength  reliability  evaluation  procesa.  However,  it  should  be 
evident  that  the  stress  analysis  procedures  are  basically  deterministic  in  nature.  Hence, 
the  variation  between  demonstrated  strength  and  calculated  strength  is  not  truly  a 
statistical  quantity  but  rather  u problem  in  analysis  procedure.  Such  miscalculations  as 
an  inadequate  assessment  of  local  buckling  or  net  section  area  of  a structural 
component  are  not  a result  of  the  material’s  structural  behavior.  Fortunately,  the 
design  and  production  of  an  aircraft  structural  system  have  usually  included  a 
demonstration  ur  verification  phase,  the  static  test,  which  effectively  ferrets  out  such 
errors,  if  existing.  Local  structural  revisions  can  correct  the  static  analysis  deficiency 
usually  in  some  reasonable  and  practicol  manner.  On  the  other  hand,  to  treat  the  entire 
structure  with  its  various  load  conditions  resulting  from  random  variations  and  usage, 
analyzing  all  possible  local  deficiencies  would  be  virtually  impossible  and  would  add  a 
structural  weight  penalty.  Furthermore,  this  weight  penalty  cannot  always  be  hilly 
compensated  with  an  increased  payload  or  performance. 

A similar  difficulty  exists  in  the  case  of  fatigue-critical  structure.  Fatigue  performance, 
whether  in  terms  of  crack  initiation,  crack  growth,  and/or  residual  strength,  is  a 
variable  resulting  from  a highly  individualistic  behavior  of  a particular  bit  of  material 
in  a specific  structural  configuration  and  variable  loading  environment.  Nevertheless, 
the  probable  behavior  of  such  material  fatigue  performance  will  respond  to  the 
nondotcrministic  approach  offered  by  reliability  analysis  methods.  Unfortunately, 
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fatigue  performance  ie  more  sensitive  to  actual  local  geometry  and  repeated  loading 
conditions  than  the  comparable  elements  of  static  strength.  Furthermore,  experience 
shows  that  deterministic  analysis  to  predict  and  interpret  fatigue  performance  is  a more 
diffloult  task,  since  it  is  much  less  understood  than  that  of  static  strength  design. 
Nevertheless,  it  is  suggested  that  analysis  procedures  for  fatigue  assessment 
incorporate  some  deterministic  aspects  as  well  as  statistical  ones  despite  the 
significauce  of  unexpected  incidents  which  occur  in  new  and/or  different  structural 
systems.  However,  careful  and  knowledgeable  treatment  of  the  elements  causing  fatigue 
damage  can  aid  in  the  precision  of  the  deterministic  phase  of  the  task.  Then  further 
assurance  is  obtained  by  recognizing  the  elements  of  variability  and  characterizing 
their  behavior  so  that  the  likelihood  of  damage  occurrence  is  a remote  or  acceptable 
possibility. 

Inherent  in  the  background  of  this  development  is  the  importance  of  structural 
configuration  and  material  selection  to  fail-safe  (i.e.,  damage-tolerant)  concepts  and  the 
associated  inspection  maintenance  plans.  This  consideration  provides  security  in 
maintaining  structural  safety  and  reduces  fatigue  damage  to  an  economic  consideration 
rather  than  one  of  safety. 


In  continuation,  the  important  feature  in  the  achievement  of  staiio  structural  integrity 
in  aircraft  structures  has  been  the  practice  of  design  and  test  verification  before  a fleet 
has  advanced  significantly  through  a production  line.  Fatigue  performance  is  responsive 
to  similar  treatment  but  with  some  reservations.  The  well-practiced  and  demonstrated 
procedures  of  static  stength  design  are  not  sensitive  to  the  local  strain  behavior  of 
actual  complex  structures.  Likewise,  the  expected  use  may  be  established  for  design 
purposes,  but  both  the  actual  use  and  response  of  the  structural  system  may  depart 
i>om  the  design  presumptions.  However,  the  analysis  difficulties  can  bo  essentially 
resolved  by  specific  developmental  and/or  verification  testing  or  by  the  use  of  previous 
experience.  Such  representative  testing  provides  an  estimate  of  fatigue  life  distribution 
location  parameters,  while  reasonable  interpretation  of  other  comparable  data  provides 
guidance  for  estimating  the  necessary  shape  parameters. 

The  condition  of  u structural  system  can  only  bo  determined  by  inspection.  The  success 
of  the  inspection  plan  depends  upon  its  sensitivity  and  thoroughness.  The  material 
properties  of  crack  growth  rate  and  residual  strength  are  the  moot  important 
parameters  in  the  formulation  of  any  such  surveillance  scheme.  The  intervals  between 
the  time  at  which  a crack  is  detectable  and  the  time  at  which  the  crack  has  developed  to 
the  critical  fail-safe  size  provides  the  only  opportunity  to  detect  the  damage.  The 
principles  of  iVacture  mechanics  provide  a means  to  estimate  crack  growth  and  resultant 
strength  of  the  damaged  structure.  Furthermore,  the  generation  of  crack  growth 
information  is  a far  simpler  task  than  determining  time  to  initiation  of  a detectable 
crack  in  structural  details.  In  addition,  the  rapid  development  of  fracture  mechanics 
information  has  provided  usable  data  on  crack  growth  rate  and  residual  strength. 
Hence,  appropriate  models  can  be  selected  and  subsequently  verified  experimentally.  As 
noted  in  earlier  studies,  the  intercept  and  slope  both  appear  to  vary  when  plotted  as  a 
function  of  crack  growth  rate  versus  stress  intensity  iactor  within  the  usef^ul  range  of 
stress  intensity  values. 
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In  aummary,  the  key  materiala/structurea  fatigue  performance  input  parameters  in  this 
proposed  application  of  reliability  methodology  include  the  following: 

1.  Time  to  detectable  crack  initiation  with  its  variability  identified  in  a distribution 
with  its  parameters  detennined  by  the  specific  detail 

2.  Crack  growth  rate  with  its  variability  identified  by  a distribution  with  its 
parameters  detemined  by  the  specific  material  of  construction 

3.  Residual  strength  or  fracture  toughness  properties  of  the  detail 

Closely  associated  with  these  material/structuros  variability  behaviors  are 
complementary  deterministic  analyses  or  supportive  programs; 

1.  Cumulative  damage  analyses  for  determining  crack  initiation  time  and  its 
responsiveness  to  environmental  exposure  of  the  particular  detail 

2,  Cumulative  crack  growth  analyses  for  the  environmental  exposure  of  the 
particular  detail 

These  input  functions  or  parameters  are  referenced  in  terms  of  specific  details.  These 
details  may  be  a particular  fitting,  a joint,  a component,  or  a typical  structural  area 
such  as  a mechanically  fastened  skin-atiffenor  structure  or  a bonded  structure. 
However,  all  the  concepts  are  developed  with  consideration  of  the  properties  of  metal 
alloy  structures  having  local  fatigue-critical  details. 


8.0  CONCLUSIONS  AND  RECOMMENDATIONS 


B.1  CONCLUSIONS 

There  are  several  goalu  of  the  statistical  investigation  that  have  been  acoompliahocl: 

1.  The  study  of  a symmetric  class  of  distributions  for  application  as  a model  for  the 
variability  of  fatigue  damage  initiation 

2.  The  development  of  unbiased  estimation  procedures  for  the  parameters  of  the 
distribution  using  censored  date 

3.  The  development  of  a technique  for  testing  the  hypothesis  of  two  censored  samples 
belonging  to  a common  parent  population 

' \ 

Additionally,  developmental  work  pertaining  to  the  reliability  of  airplane  structures 
with  mixed  populations  resulting  from  different  lovels  of  fatigue  performance  capability 
has  been  incorporated  into  the  structural  reliability  analysis  system.  Several 
parametric  evaluations  of  diverse  structural  variables  were  conducted  using  the 
reliability  analysis  system,  and  the  results  have  been  presented  and  discussed. 

Conclusions  that  may  be  drawn  fVom  this  investigation  include  the  following; 

1.  The  symmetric  throe-parameter  distribution  does  describe  the  empirically 
determined  fatigue  distribution. 

2.  I'he  flexure  parameter  of  this  distribution  may  be  presumed  to  be  a constant,. 

3.  MLE’s  of  the  remaining  parameters  are  computable  from  censored  data. 

4.  A practical  simulation  procedure  for  censored  data  has  been  duvolopud,  an 
Important  step  resulting  in  the  capability  of: 

e Computing  bounds  for  censored  data 

• Obtaining  unbiased  point  estimates  from  censored  data 

• Testing  the  hypotheses  that  two  independent  censored  samples  have  a 
common  parent  population 

5.  Structural  reliability  is  BlgnificanUy  affected  by  the  magnitude  and  the  tVequenoy 
of  initial  flaws  or  damage. 

Parametric  studies  conducted  in  this  investigation  have  shown  that  structural 
reliabiltiy  is  dependent  on  several  interdependent  design  variables,  such  as; 


• Loading  environment 

• Residual  strength 

• Design  configuration 

• Material  characteristics 

Finally,  the  reliability  model  with  the  concept  of  initial  flaws  or  mixed  populations  of 
fatigue  performance  can  reflect  the  extreme  initial  behavior  shown  in  the  collected 
groups  of  experimental  fatigue  test  data.  Therein,  the  empirical  cumulative  failure 
distribution  when  plotted  on  normal  probability  paper  exhibits  a characteristic  flexure, 
If  there  is  a positive  probability  of  an  initial  crack  at  time  aero,  the  cumulative 
distribution  will  become  asymptotic  to  this  value.  Furthermore,  if  there  is  a waiting 
time  until  the  crack  initiates,  the  cumulative  distribution  function  will  rise  more 
steeply  in  the  early  extreme  portion  than  the  empiric  fatigue  performance  distribution 
function.  Hence,  inspection  and  maintenance  control  plans  should  be  based  on  models 
that  reflect  such  conditions,  and  not  upon  linear  extrapolations  from  data. 


8.2  RECOMMENDATIONS 

It  is  recommended  that: 
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1.  The  reliability  analysis  system  be  expanded  to  consider  more  complex  structural 
configurations  and  loading  conditions.  An  Initial  attempt  to  evaluate  a 
mechanically  fastened  skin-stiffener  type  structure  is  described  in  reference  14, 
and  some  results  are  presented.  It  is  proposed  that  this  technology  be  developed  to 
encompass  flaws  and  fatigue  cracks  in  either  the  skin  or  the  stiffener,  at  or  away 
fVom  the  fastener  holes,  under  typical  flight  loads. 

2.  A corrosion  term  bo  included  in  the  reliability  analysis  system  to  complement  the 
initial-flaw  model  as  u possible  source  of  early  fatigue  incidents. 

3.  A capability  to  modify  selected  input  parameters  during  the  operational  lifetime  of 
a structure  be  incorporated  into  the  reliability  analysis  systtm.  This  would 
facilitate  the  application  of  the  analysis  system  to  fleet  evaluation  and 
management  tasks. 

4.  The  collection  and  evaluation  of  suitable  infomiation  be  undertaken  to  determine 
the  magnitude  and  fVeqviency  of  such  preflaws  typical  of  current  airplane 
structures. 

6.  The  data  simulation  procedures  be  used  to  develop  a usable  matrix  of  factors  which 
can  be  used  to  unbias  estimates  from  censored  samples. 

6.  The  statistical  methods  developed  for  application  to  censored  data  be  extended  to 
include  other  distribution  models. 
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Table  3.~-Comparison  of  Parameters  From  Data  With  Differing  Life  Lengths 


Table  6.— Scale  and  Location  Parameter  Estimates  of  Censored  Test  Data  Samples 


Sample  no. 


MLE 

scale  parameter 

MLE 

location  parameter 

Scale  parameter  ratio 
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12.692 
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Table  7.— Location  Parameter  Estimates  of  Paired  Censored  Test  Data  Samples 
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Comparison  of  Aiuminum  Ailoy  Fatigue  Test  Data  Distribution  and  the 
Birmbaum-Saunders  Inverse  Hyperbolic  Sine  Distribution 
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Figure  7.— Iteration  Resu/ti  of  Fatigue  Data  Sanyjle,  Open  Hole  Specimens  Randomly 
Loaded— ^mph  Size  = 10 
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Figure  8. -Comparison  of  Aluminum  Alloy  Fatigue  Test  Specimen  Types 


Cumulative  -^equency  Cumulative 


12.54 


tZ44 


t-0 


ljr-° 


Ceniored  fatigue  data  sample 


16  failure 
observations 

32  suspended 
observations 


132  000 

143  000 

147  000 

168  000 

150  000 

162  000 

181  000 

191  000 

192  000 

204  000 

236  000 

240  000 

241  000 

244  000 

246  000 

248  000 

260  000 
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Ftgur0  12,~Simu(0tlon  of  Cansot»d  Sample  With:  4 1 Failure  Observations;  1 1 Censored  Observations 
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Figure  13.-Slmulat(on  of  Censored  Sample  With;  33  Failure  Observations;  20  Censored  Observations 
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Figure  14,—Simulaticii  of  Cfiisored  Sample  With;  8 Failure  Observations;  40  Censored  Observations 
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'Simulation  of  Censored  Samples  With  Various  Assumed  Flexure  Parameters 


rgt/re  18.— Example  Curve  of  Unbiasing  Factors  for  the  Location  Parameter 
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Figure  20- Distribution  of  Parameter  Statistics  for  Two  independent  Censored  Data  Sainpies- 
Sample  1:  33  Failures,  20  Suspensions;  Sample  2:  33  Failures,  20  Suspensions 
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Figure  21.— Distribution  of  Parametor  Statistics  for  Two  indopendent  Censored  Data  Samples 
Semple  1:  10  Failures,  44  Suspensions;  Sample  2:  29  Failures  2B  Suspensions 
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These  curves  ore  based  on  data  recorded  from 
approximately  18  000  hr  of  operations  of 
four-engIned  jet  transport  airplanes 


Frequency  exceedances/hr 


FIgun  22.— Exceedances  of  Gust-Induced  Load  Ratios  for  Two  Critical  Major  Locations  of  a 
Jet-Engined  I anker /Transport-Type  Airplane 
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Figure  24.-1  mpact  of  Residual  Strength  on  the  Reliability  of  a Skin-Type  Structure 
Without  Initial  Flam-Type  A Loads 
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Figure  2S.~lmpect  of  Residua!  Strength  on  the  Reliabliity  of  a Skin-Type  Structure 
With  initiai  Flaws—  Type  A Loads 
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